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ABSTRACT 


Experimental  damage  thresholds  for  Si,  Ge , GaAs,  GaAs^P^^ 
and  GaP  have  been  determined  at  10.6  pm  and  1.06  ym  laser  wave- 
lengths. For  CW  CC>2  laser  damage,  all  the  thresholds  are  of  the 
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order  of  a few  kW/cm  for  an  irradiation  time  of  about  0.10  second 

At  1.06  ym  wavelength  and  for  nanosecond  pulses,  the  damage  thresh 
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olds  are  of  the  order  of  10  W/cm  . The  damage  thresholds  can  be 
explained  by  assuming  that  the  anomalous  absorption  takes  place 
through  the  parametric  instability  of  the  type  formulated  by 
DuBois  and  Goldman.  The  plasma  concentration  required  for  the 
onset  of  this  instability  at  10.6  ym  is  easily  achieved  by  photon 
induced  excitation  of  the  valence  electrons  to  the  conduction  band 
through  a continuous  distribution  of  Shockley  surface  states.  At 
1.06  ym,  the  required  plasma  concentration  in  case  of  Germanium 
and  Silicon  is  produced  by  band  to  band  transitions  and  by  Zener 
tunneling  in  case  of  GaAs,  GaAs^P.^^  and  GaP. 

In  practice,  the  damage  occurs  both  at  the  surface  of  the 
semiconductor  and  within  the  bulk,  creating  defects  that  cause 
carrier  removal  and  mobility  degradation.  The  defects  so  created 
are  more  active  at  lower  temperatures  than  at  room  temperature 
as  seen  from  the  experimentally  measured  Hall  mobility  versus 
temperature  curves  before  and  after  laser  damage.  Carrier  re- 
moval and  mobility  degradation  curves  appear  qualitatively  to 
be  similar  to  those  observed  for  electron,  neutron  and  y-ray 
damage.  These  changes  in  electrical  properties  can  be  explained 
by  assuming  the  creation  of  localized  traps  in  the  band  gap  and 
applying  the  model  of  James  and  Lark-Horovitz . 
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CHAPTER  1 
INTRODUCTION 

The  study  of  the  numerous  ways  in  which  radiation  and  matter 
can  interact  has  been  a subject  of  intensive  research  for  several 
years.  The  interaction  can  be  linear  or  nonlinear  depending  upon 
the  properties  of  matter  and  radiation.  The  end  result  of  the 
interaction  is  mechanical  or  electrical  damage  to  the  matter. 

There  have  been  numerous  studies  of  radiation  induced  damage 
in  semiconductors  due  to  electrons  [1,  2,  3],  neutrons  [4,  5,  6], 
deutrons  [7,  8,  9],  a-particles  [10,  11]  and  y-rays  [2,  12].  The 
study  of  radiation  induced  damage  in  electrical  properties  of  semi- 
conductors is  of  great  importance  and  is,  therefore,  a subject  of 
continuous  research  in  solid  state  electronics. 

The  present  study  of  laser  induced  damage  in  semiconductors 
was  undertaken  for  several  important  reasons.  Unlike  the  extensive 
data  available  on  particle  and  y-ray  damage,  a literature  survey 
shows  that  only  meager  and  incomplete  data  exist  on  laser  damage 
of  semiconductors.  Today  more  and  more  semiconductor  devices  are 
finding  applications  in  a wide  variety  of  fields  such  as  in  laser 
systems,  laser  energy  conversion  [13],  air  force  weapons,  satel- 
lites (exposed  solar  cells),  IR  controlled  weapons,  opto- 
electronics and  micro-optics.  Therefore,  a knowledge  of  the  be- 
havior of  the  electronic  properties  of  semiconductors  under  the 
effect  of  laser  radiation  is  essential. 


* 

Numbers  in  brackets  refer  to  references  at  the  end  of  each 
chapter. 
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Next,  the  great  structural  sensitivity  of  semiconductors  to 
radiation  makes  them  very  useful  systems  in  which  to  investigate 
the  nature  of  laser  damage  in  solids.  The  observed  laser  induced 
changes  in  electrical  properties  such  as  carrier  mobility,  life- 
time and  concentration  can  be  used  to  diagnose  the  damage  mechan- 
isms that  cause  catastrophic  material  failure  in  optical  materials 
for  high  power  lasers.  For  obvious  reasons,  such  techniques  can- 
not be  applied  to  metals  and  insulating  dielectrics. 

The  damage  to  semiconductor  materials  can  occur  through  a 
number  of  mechanisms.  Chapter  2 describes  some  of  the  more  im- 
portant mechanisms  through  which  a laser  beam  can  interact  with 
the  semiconductor  matter.  The  absorption  of  energy  can  take 
place  through  linear  processes  like  band  to  band  transitions, 
free  and  bound  carrier  excitation,  lattice  absorption  or  it  can 
occur  through  nonlinear  processes,  e.g.,  avalanche  mechanism, 
self-focusing  or  plasma  instabilities.  Most  of  these  mechanisms 
are  described  in  detail  in  this  chapter. 

Chapter  3 gives  the  theory  of  the  Hall  effect  as  applied 
by  Van  der  Pauw  to  determine  the  carrier  mobility  in  semicon- 
ductors. The  experimental  details  regarding  the  fabrication  of 
devices  and  the  experimental  setup  to  determine  mobility,  car- 
rier concentration  and  conductivity  are  given.  Since  radiation 
damage  can  alter  the  lattice  structure  and  introduce  extra- 
scattering centers,  the  theory  of  different  types  of  scattering 
mechanisms  is  also  included  in  this  chapter. 

Experimental  results  on  CW  CC>2  laser  damage  in  different 


2 


semiconductors  are  given  in  Chapter  4.  It  is  shown  that  the  ex- 
perimental damage  thresholds  can  be  explained  [14]  by  applying 
the  theory  of  parametric  instability  as  formulated  by  DuBois  and 
Goldman  [15] . The  observed  mobility  and  carrier  concentration 
changes  are  explained  [16]  by  applying  the  model  of  James  and 
Lark-Horovitz  [17] . 

Chapter  5 gives  the  experimental  results  on  damage  thresholds 
for  semiconductors  for  Nd : YAG  laser  (wavelength  1.06  x 10  m) . 

The  theory  of  parametric  instability  again  explains  the  damage 
thresholds  quite  well.  It  is  shown  that  the  results  of  modified 
two  beam  instability  of  Kaw  and  Dawson  [18]  do  not  explain  the  ob- 
served damage  thresholds.  The  theory  of  photon -pi asmon  insta- 
bility of  Gersten  and  Tzoar  [19]  is  also  applied  to  calculate  the 
damage  thresholds  for  different  semiconductors,  and  it  is  seen 
that  the  theory  does  not  predict  the  damage'  thresholds  correctly. 
These  two  instabilities  have  been  used  incorrectly  by  Shatas  et  al . 
[20]  to  explain  the  damage  thresholds  for  GaAs  at  CO^  laser  wave 
length. 
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CHAPTER  2 

FUNDAMENTAL  DAMAGE  MECHANISMS 

This  chapter  describes  the  possible  mechanisms  via  which 
electro-magnetic  radiation  energy  can  be  absorbed  by  matter,  par- 
ticularly by  semiconductors.  When  the  energy  absorbed  by  a semi- 
conductor exceeds  a certain  level,  mechanical  or  electrical  damage 
can  occur.  This  energy  is  called  the  damage  threshold  of  that 
semiconductor.  Here  mechanical  damage  means  any  sign  of  burning 
of  the  surface  or  pit  formation  or  any  other  visible  change  in 
the  appearance  of  the  semiconductor  surface.  Electrical  damage 
stands  for  any  changes  in  carrier  mobility,  carrier  concentra- 
tion or  carrier  lifetime. 

The  following  damage  mechanisms  will  be  described: 

1.  Band  to  band  transitions 

2.  Free  carrier,  bound  carrier  and  lattice  absorption 

3.  Avalanche  mechanism 

4.  Self  focusing 

5.  Plasma  instabilities. 

Here  it  should  be  pointed  out  that  this  list  is  by  no  means 
intended  to  be  complete  and  exhaustive. 

1.  BAND  TO  BAND  TRANSITIONS  [1] 

The  probabilities  associated  with  band  to  band  transitions 
were  first  calculated  by  Kane  [2]  in  1962,  and  the  theoretical 
predictions  were  promptly  verified  by  Allen  and  Gobeli  [3].  The 
absorption  through  a transition  from  valence  band  to  conduction 
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band  of  a semiconductor  is  called  fundamental  absorption. 


Two  types  of  transitions  must  be  distinguished.  Those  in 
which  only  a photon  is  involved  are  called  direction  transitions, 
while  those  in  which  both  a photon  and  a phonon  are  involved  are 
called  indirect  transitions.  In  indirect  transitions,  phonons 
help  conserve  the  momentum.  Below,  direct  transitions  in  both 
direct  and  indirect  band  gap  semiconductors  will  be  considered. 

Direct  Transitions  (Direct  Band  Gap  k = k . ) 

• max  mm 

In  any  transition  crystal  momentum  must  be  conserved.  Hence 
in  a direct  band  gap  semiconductor  only  vertical  transitions  will 
be  allowed  (Figure  2-1) . The  minimum  frequency  of  the  photon  that 


J 


can  induce  such  a transition  is  given  by 


hv  = E (k)  - E (k)  (2.1) 

c v 

Assuming  that  the  conduction  band  stays  flat  near  k = 0, 

then  E (k)  in  the  vicinity  of  k = 0 can  be  taken  as  a constant, 
c 

Differentiating  Equation  2.1,  we  get 

hdv  = - dE  (k) 
v 

= - d(-E  - E') 

g 

or  hdv  = dE ' (2.2) 


duction  band  and  E'  from  the  top  of  the  valence  band. 


Now  the  number  of  transitions  Nfcdv  in  the  interval  dv  can  be 
written  as 

N dv  = P N (E 1 ) dE  ' (2.3) 

t k 

where  is  the  probability  associated  with  the  transition  and 

N(E')dE'  is  the  density  of  states  in  the  energy  interval  E'  and 
E'  + dE ' . For  a parabolic  band 

N ( E ' ) = a/e7  (2.4) 

where  A is  a constant  determined  by  the  nature  of  the  semi- 
conductor . 


From  Equations  2.3  and  2.4,  one  can  write 


stfn-  - 


a,  = B(hv  - E ) **  hv  > E 
d g g 


} (2.7) 

= 0 hv  < E \ 

g ' 


The  value  of  B has  been  found  quantum-mechanically  by  Bardeen, 


Blatt  and  Hall  [4]  and  is  given  by 


ire2  (2m  ) 3//2 

B = — f (2.8) 

nch2me 


where  n is  the  real  part  of  the  refractive  index  and  mf  is  the 
reduced  mass  of  the  electron  effective  mass  and  hole  effective 


mass  given  as 


f.  _ is  called  oscillator  strength  and  is  usually  of  the  order  of 
if 

unity . 

If  m*  * m*  = m,  then  for  n = 4 
e h 

a , = 6.7  x 104(hv  - E ) **  f.,(cm-1).  (2.9) 

d gif 

Direct  Transitions  (Indirect  Band  Gap  k / 

max min 

The  band  structure  of  an  indirect  band  gap  semiconductoi 
is  shown  in  Figure  2-2. 


Valence  Band 


Figure  2-2.  A direct  (A1)  and  an  indirect  (B1) 
transition  in  an  indirect  band 
gap  semiconductor 

In  this  case  the  maxima  of  the  valence  band  and  minima  of 
the  conduction  band  lie  at  different  values  of  k. 

The  absorption  curve  for  Ge  has  been  obtained  by  Macfarlane 
et  al.  [5]  and  is  shown  in  Figure  2-3.  The  peak  in  the  absorp- 
tion curve  is  due  to  a vertical  transition  of  the  type  A'  and  an 
exciton  line.  The  absorption  at  longer  wave  lengths  is  due  to 
the  indirect  transitions  of  the  type  B'  and  the  momentum  is  con- 


served by  emission  or  absorption  of  one  or  more  phonons. 


Figure  2-3.  Absorption  in  Ge  due  to  vertical 
transitions  at  k = 0 


Indirect  Transitions  (k  . k ) 

min max 

As  noted  earlier,  in  indirect  transitions  one  or  more  phonons 

must  be  involved  to  conserve  the  crystal  momentum.  Therefore,  the 

minimum  photon  energy  required  for  such  transitions  is  given  by 

hv  = E - E (2.10) 

g P 

or  hv  = Eg  + Ep  (2.11) 

X 

Equation  2.10  implies  absorption  while  Equation  2.11  implies 

emission  of  a phonon  in  the  process. 

•n 
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The  absorption  coefficients  for  such  transitions  have  been 
theoretically  obtained  by  Bardeen  et  al.  [5] . Though  their  for- 
mula does  not  separate  the  contributions  arising  from  the  emis- 
sion or  absorption  of  phonons,  it  does  take  into  account  the  two 
possible  ways  in  which  a transition  from  a state  in  the  valence 
band  to  a state  in  the  conduction  band  can  take  place.  This  is 
shown  in  Figure  2-4. 


Figure  2-4.  Indirect  transitions  from  the  valence 
band  to  the  conduction  band 

Transition  1 corresponds  to  an  electron  excited  from  the 
valence  band  to  the  conduction  band  without  significant  change 
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in  crystal  momentum.  Since  in  this  state  of  the  conduction  band 


the  electron  has  higher  energy,  it  makes  a transition  2 to  the 
state  of  lower  energy  in  the  conduction  band  with  the  emission 
of  a phonon  of  momentum  k^  . Alternatively  an  electron  may  be 
excited  from  a state  in  a valence  band  to  a state  in  the  con- 
duction band  through  transition  3.  The  hole  left  in  the  valence 
band  then  makes  a transition  to  a state  near  the  maxima  of  the 
valence  band  by  emission  or  absorption  of  a phonon. 

To  determine  the  total  absorption  coefficient  one  must  add 
the  contribution  from  the  transitions  involving  the  absorption 
of  a phonon  or  emission  of  a phonon. 


aT  = a + a, 
I E A 


(2.12) 


where  is  the  total  absorption  coefficient  due  to  indirect 

transitions,  a and  a,  are  the  absorption  coefficients  due  to  in- 

E A 

direct  transitions  involving  the  emission  and  absorption  of  a 
phonon  respectively. 

From  Equations  2.10  and  2.11,  we  see  that 

ci  = 0 for  hv  < E + E (2.13) 

E 9 P 

a =0  for  hv  < E - E (2.14) 

A g P 

If  k = 0 and  k . are  well  separated,  then  the  transition 
mm 

probability  between  a state  near  the  top  of  the  valence  band  and 
a state  near  the  bottom  of  the  conduction  band  will  vary  only 
slightly  with  k and  may  be  considered  as  a constant.  The  total 
absorption  coefficient  will  then  be  determined  by  the  density  of 
states  from  which  the  transitions  are  taking  place  and  the  relative 


1 3 


probability  of  emission  or  absorption  of  a phonon  of  energy  E . 


N 


The  energy  E will  be  nearly  Rk  . and  can  be  treated  as  constant, 
p min 

The  conservation  of  energy  for  such  transitions  requires 

hv  = E ± E + E + E ' (2.15) 

g P 

where  E is  the  energy  of  the  electron  above  the  conduction  band 
and  E'  is  the  energy  below  the  valence  band.  If  E'  is  fixed, 
then 

hdv  = dE  (2.16) 

The  number  of  states  in  the  conduction  band  in  this  energy 
interval  is  given  by 

N (E)dE  = aE  dE 
c 

= a (hv  - E ± E - E'^dE  (2.17) 

g P 

where  a is  a constant. 

Now  the  total  number  of  the  pairs  of  states  between  which 

the  transitions  lying  in  the  interval  v and  v + dv  are  allowed 

can  be  found  by  integration  over  all  E'  satisfying  Equation  2.15. 

From  Equation  2.15,  the  maximum  value  of  E'  is  given  by 

E'  = hv  - E + E . (2.18) 

m g P 

For  the  density  of  states  of  the  valence  band  in  the  inter- 
val dE'  we  have  a'E'^dE'.  Hence  the  required  number  of  pairs  of 
states  is  given  by  integrating  Equation  2.17 

E* 

^ L L 

N ( v) dv  = aa'hdv  f (E'  - E')  E’  dE* 

m 
o 

= DE' 2dv  (2.19) 

m 

where  D is  a constant. 
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For  the  absorption  of  a phonon 


E'  = hv  - E + E 

m q p 


and,  therefore. 


N(v)dv  = D(hv  - E + E ) dv 
g P 


(2.20) 


Phonons  are  Bose  particles  and  the  probability  that  a phonon 
will  have  an  energy  E^  is  given  by 


lIp^Ep^  exp(E  /k  T)  - 1 
P B 


(2.21) 


so  that  the  absorption  coefficient  aft  is  given  as 


A' (hv  - E + E ) 

a = 2 P_ 

A exp(E  /k  T)  - 1 
P B 


for  hv  > E - E 
g P 


for  hv  < E - E 

g p 


(2.22) 


where  A'  is  a constant. 


Similarly  we  can  show  that 


A' (hv  - E - E ) 

a = 2 E_ 

E 1 - expf-Ep/k^T) 


for  hv  > E + E 

g p 


for  hv  < E + E 

g P 


(2.23) 


(hv  - E - E ) (hv  - E + E ) 

a = A* + E — hv  > E + E 

I 1 - exp(-E  /k  T)  exp (E  /k  T)  - 1 g P 

p B P B 


(hv  - E + E 

= A. 5 P— 

exp(E  /k)  - 1 
p BT 


E - -E  < hv  < E + E > (2.24) 

g p g P l 


hv  < E - E 

g P 


v -w.; 


The  phonon  considered  in  the  above  treatment  can  either  be  an 


r 


a 


acoustic  or  an  optical  phonon.  In  each  mode  it  can  either  be  a 

longitudinal  or  transverse  type.  This  will  split  the  absorption 

curve  into  four  parts.  Using  the  analysis  given  above,  the  values 

of  longitudinal  and  transverse  acoustic  phonon  energies  for  Ge 

have  been  found  [6]  to  be  as  E = 0.008  eV,  and  E „ = 0.027  eV. 

pt  pd 

In  the  case  of  Ge  the  absorption  due  to  optical  phonons  is  negli- 
gible . 

For  Si,  both  acoustic  and  optical  phonons  give  appreciable 
contribution  to  the  absorption.  The  phonon  energies  have  been 
found  to  be 

Acoustic  Phonons  E . = 0.0185  eV,  E „ = 0.0575  eV 
pt  pd 

Optical  Phonons  E . = 0.120  eV,  E „ = 0.091  eV 
pt  pd 

A more  exact  treatment  of  absorption  due  to  direct  and  in- 
direct transitions  has  been  given  by  Fan  [7] . 


Transitions  Involving  Surface  States 

If  the  surface  state  density  is  high  and  the  lifetime  of 
these  states  is  long,  then  an  appreciable  absorption  can  take 
place  through  surface  states.  The  quantum  yields  for  such  pro- 
cesses have  been  calculated  by  Kane  [2]  and  can  be  very  high. 


2.  FREE  CARRIER  ABSORPTION  [18] 

When  the  incident  photon  energy  is  not  enough  to  cause  band 
to  band  transitions,  energy  can  still  be  absorbed  due  to  intra- 
band transitions.  This  is  called  free  carrier  absorption.  Ab- 
sorption due  to  intraband  transitions  is  proportional  to  the 
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number  of  free  carriers  in  the  band. 


In  the  non-degenerate  case,  semiclassical  treatment  can  be 


applied  to  the  problem  of  free  carrier  absorption.  The  equation 
of  motion  of  a free  carrier  in  an  alternating  electric  field  of 
the  incident  wave  can  be  written  as 


m*v(t)  = -eE(t) 


(2.25) 


where  m*  is  the  effective  mass  of  the  free  carrier  and  v its  ve- 
locity  under  the  action  of  the  applied  delectric  field  E. 

In  the  presence  of  randomizing  collisions  Equation  2.25  will 
be  modified  to 


(v(t)  + XilL)  = eEe" 


* (2.26) 


where  t is  the  appropriate  relaxation  time  of  the  carriers  in  the 
semiconductor.  The  charged  carrier  under  the  action  of  the  ap- 
plied harmonic  field  Ee  1Wt  will  be  forced  to  oscillate  with  a 
velocity 

->  -iut 
v ( t ) ~ ve 


Substituting  this  in  Equation  2.26  and  solving  for  v,  one 


(2.27) 


m*  (—  - ito) 
't 


The  current  density  can  now  be  written  as 


j = nev 


_ -n  e - t t 

m* ( 1 - ion) 


(2. 20) 


where  e is  the  charge  on  the  carrier  (electron  or  hole) . 


From  Equation  2 .28,  one  can  write  for  the  conductivity  as 


o (to)  = 


ne2T  / 1 \ 

m*  \ 1 - itoi  / 

o [ r \ 

oil-  110T  I 


(2.29) 


where  0q  is  the  dc  conductivity  of  the  carrier  species. 

Following  a similar  procedure  and  assuming  r - r^e  10jt/  one 
can  solve  Equation  2.26  for  r and  then  write  the  polarization  vec- 


P = ner 


2 

m*  (to  + iio/t) 


Now  the  dielectric  constant  e(io)  is  given  by 


£ (to)  = 1 + - 
E 


= 1 + 


-m*(w  + iio/x) 


£ (to)  = 1 - — 


t (to  + ito/t) 

Rewriting  Equation  2.31, 


E((0)  = i 


\ 1 + to  TZ/  to(l  + 


, 2 2 

+ (0  T 


(2.31) 


(2.32) 


The  real  part  of  this  dielectric  constant  when  equated  to 
zero,  gives  the  resonance  frequency  of  the  carrier  plasma  and  can 


be  shown 


)4ti  ne' 

y m* 


(2.3  I) 
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The  imaginary  part  gives  the  damping  rate  of  the  electromag- 


netic wave  or  the  absorption  coefficient,  which  in  this  case  is 

1 


o i 
o 


, 2 2 I nee 

I+wt/  O 


2.34) 


where  n is  the  real  part  of  the  refractive  index,  c the  velocity 
of  light  and  eq  is  the  permittivity  of  free  space. 

Now  two  cases  arise: 

(a)  ojt  > > 1,  which  is  true  for  most  IR  frequencies 
o 


\ 0)  T / O 


Substituting 

m*u 


and 


o = neu 
o c 


a = 


2 *2 

oj  m*  p nee 
c o 


(2.35) 


This  equation  shows  that  the  absorption  coefficient  for  free 
carriers  is  directly  proportional  to  the  carrier  concentration 
and  the  square  of  the  wavelehg-th. 

(b)  For  very  low  frequencies  o)T  < < 1 
and,  therefore 


ncr. 


3.75  (~j  0o  cm  ^ 


(2.  ?<) 


where  l is  the  dielectric  constant  of  the  semiconductor. 
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Expressions  2.35  and  2.36  have  been  experimentally  verified  by 
by  Gibson  [9) . 

Bound  Carrier  Excitation 

Materials  that  are  characterized  by  a band  gap  show  a low 
frequency  tail  of  the  fundamental  absorption.  This  absorption  co- 
efficient can  be  expressed  as 


Figure  2-5.  Optical  absorption  versus  optical  frequency 

for  a typical  material  (After  Reference  [10]) 
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(2.37) 


a (i*J)  = A exp  [B  (hid  - C )/k  T] 

9 tJ 

where  A and  £3  are  constants.  This  expression  is  called  Urbach's 
lule  [10,  11).  There  have  been  several  theoretical  attempts  to 
explain  this  result  and  none  of  them  is  satisfactory.  The  most 
recent  attempt  to  explain  this  effect  is  by  Dow  and  Redfield  [12 J. 

The  Urbach  tail  is  usually  very  steep  and  can  produce  signifi 
cant  absorption  for  those  materials  whose  band  gap  is  not  more 
than  20  k T above  the  frequency  which  one  wants  to  transmit  [13] . 
This  is  shown  in  Figure  2-5. 

Lattice  Absorption 

The  vibration  spectrum  of  diatomic  lattices  has  two  branches, 
namely  acoustic  and  optical.  The  electric  field  of  the  incident 
electromagnetic  wave  can  excite  optical  vibrations  in  the  lat- 
tice, thereby  causing  strong  absorption  at  the  resonance  frequency 
A number  of  investigations  have  shown  that  the  absorption  coef- 
ficient can  be  written  as  [14] 

-w/u> 

a (id)  = Ae  ° (2.3U) 

where  is  of  the  order  of  Restrahl  frequency  and  varies  with 
lattice  temperature.  The  presence  of  impurities  can  cause  extra 
absorption  by  providing  different  modes  of  vibration  in  the  crys- 
tal. An  excellent  review  of  the  current  work  is  given  in 
Reference  [15]. 

3.  AVALANCHF  MECHANISM  [16,  17] 

The  approximate  equality  between  the  rms  field  of  the  l.iset 
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pulse  and  the  dc  breakdown  field  of  several  alkali  halides  and 
other  band  gap  materials  has  led  Bloembergen  [17]  to  believe  that 
the  process  of  laser  induced  breakdown  is  electron  avalanche  sim- 
ilar to  the  dc  avalanche  breakdown  of  reverse  biased  p-n  junction.-.. 

This  theory  is  an  extension  of  Shockley's  dc  avalanche  mech- 
anism. A few  "lucky  electrons"  undergo  momentum  reversing  col- 
lisions just  in  time  to  be  further  accelerated  by  the  reversed 
field  of  the  laser  pulse.  In  other  words  the  collision  time  of 
the  electron  is  an  odd  integral  multiple  of  the  half  period  of  the 
laser  pulse. 

The  equation  of  motion  of  an  electron  in  the  presence  of 
collisions  can  be  written  as 


m*r  + eEe  la>t  + v m*r  = 0 (2.39) 

e c e 

where  Ee  la)t  is  the  laser  field,  v is  the  collision  frequency  of 

c 

the  electron  in  the  crystal  lattice  and  m*  is  the  effective  mass 
of  the  electron. 

Solving  Equation  2.39  for  r produces 


• e c ( exb 

m*  ‘ \-iu>  + v 
e c 

The  rate  at  which  an  electron  moving  with  the  above 

tracts  energy  from  the  electron  is  given  by 


(2.4' 

velocity  ex- 


• 1 

W = - Re ( r 


eE*) 


m* (-iw  + v) 
e 


„ „ , 2 2 
2m*  (v  + o)  ) 
e c 


(2.41) 
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Now  writing 


v = 
c 


m*y 

e e 


therefore , 


W = 


yeE 
2 2 

(1  + u>  t ) 


(2.42) 


where  E can  be  considered  as  rms  electric  field  of  the  laser. 

When  the  energy  gained  by  the  electron  exceeds  the  band  gap 

of  the  material,  it  will  cause  further  ionization.  This  process 

will  keep  multiplying  till  the  electron  concentration  is 
17  18  ~3 

10  ~ 10  cm  , which  gives  an  exponential  contribution  to  the 


index  of  refraction 


(h  + i<) 


2 

2ne  t (i  - lot ) 

2 2 

g m*w  (1  + a)  T ) 
o e 


n exp 


t 

/ a (E) dt 
o 


(2.43) 


where  n is  the  linear  index  of  refraction  and  a (E)  is  the  elec- 


o 

trie  field  dependent  ionization  rate.  n is  the  original  number 
of  electrons  per  unit  volume.  The  exponential  absorption  even- 
tually causes  damage. 

The  application  of  this  process  to  explain  laser  damage  thr<-  ' 
olds  for  semiconductors  has  recently  been  questioned  by 
Shatas  et  al . ( 18 ] . 


4.  SELF-FOCUSING  [19,  20,  21,  22,  23] 

For  most  purposes,  the  refractive  index  of  a substance  can 
be  specified  as  a single  number.  However,  it  is  well  known  that 
the  refractive  index  is  dependent  on  temperature,  wavelength  and 
electric  field.  When  a medium  is  sufficiently  absorbing,  the 
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refractive  index  has  to  be  represented  as  a complex  number. 

An  applied  electric  field  polarizes  a medium.  The  simplest 
case  is  when  the  polarization  is  proportional  to  the  applied  elec- 
tric field.  It  can  be  mathematically  stated  as 


(2.44) 


where  x is  known  as  the  susceptibility  of  the  medium.  In  an  aniso- 
tropic medium  x may  be  different  in  different  directions  and  has 
to  be  expressed  as  a tensor  quantity.  In  Equation  2.44,  x is  in- 
dependent of  the  electric  field.  For  very  high  electric  fields, 
local  polarization  cannot  be  neglected  and,  therefore,  x would 
become  a function  of  the  applied  electric  field.  However,  the 
dependence  of  ? on  E must  be  such  that  it  reduces  to  Equation 

2.53  for  low  electric  fields.  Also  as  E reverses,  P should  re- 
_ “>• 

verse  too.  In  other  words  P can  be  a function  of  only  odd  powers 

—y 

of  E.  Hence  one  can  write, 


->nf  2 4 \ -*■ 

P = Co(Xl  + *2E  + *3E )E 


(2.45) 


Neglecting  X3E  and  higher  order  terms  in  the  nonlinear 


polarization , 


Xeff  X1  + X2E 


(2.40) 


Now  the  refractive  index  is  given  by 


neff  (1  + xef f ’ 


* (1  * x/  * \ 


(1  + Xx> 


°r  neff  = r,o  + °2E 


(2.47) 


’W.Z 


(2.48) 


0 , however,  is  limited  due  to  diffraction  and  for  a beam  of 
diameter  2k  it  is  given  as 

1.221  (2.4'<) 

d 2ns 

o 

Depending  on  the  value  of  0 , the  following  three  cases  arise. 

(a)  < 0 , 

o d 

In  this  case  the  beam  would  spread. 

(b)  6=0, 

o d 

In  this  case  the  beam  would  not  spread  and  would  get  self- 
trapped. 

(c)  0 > 0^, 

o d 

The  beam  would  converge  or  be  self  focused. 

Case  (b)  determines  the  threshold  power  for  self  focusing. 


Hence  the  threshold  power  for  self  focusing  can  be  written 


th  8 ii 

c 

64n 

For  CS2, 


- (1.22A)2 

2 


2 


ESU 


( 2.r>l ) 


then  at 

-4 

X = 0.7  x 10  cm 
P1*  = 50  kW 

Such  powers  are  easily  obtainable  in  modern  day  lasers. 
The  nonlinear  refractive  index  n2  can  be  due  to  several 
causes.  The  three  main  contributions  to  n2  come  from 

(a)  Thermal 

(b)  Electrostriction 

(c)  Kerr  Effect. 


(a)  Thermal  [24] 

Heating  of  a material  causes  density  changes  which  according 
to  Clausius  Mossoti  relation  would  change  the  refractive  index. 
Thermal  contribution  to  r\^  is  negative  and,  therefore,  causes 
defocusing  action. 

(b)  Electrostriction  [25] 

Electrostriction  is  the  mechanical  deformation  caused  by  the 
electric  field  of  the  optical  beam.  This  causes  an  increase  in 
density  and,  therefore,  an  increase  in  the  refractive  index.  Tin 
contribution  to  n2  in  this  case  is  positive  and  causes  self- 
focusing.  This  effect  is  very  important  in  solids. 
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c)  Kerr  Effect  (26] 


ft 

K 


l: 


The  increase  in  n2  due  to  Kerr  effect  is  very  important  in 


liquids  with  anisotropic  molecules.  The  optical  field  of  the 
laser  beam  produces  an  alignment  effect  on  these  molecules.  It 
can  be  shown  that 


- 27T(Pm~  Pi)No 
n2  35k  Tn 


(2.52) 


B 'o 


where 


L = 


,2  4 

(n  + 2) 
o 


81 


is  the  Lorentz  local  field  correction  factor.  N is  the  number 

o 


density  of  the  molecules,  pt  ( and  p^  are  respectively  the  princi- 
pal polarizabilities  of  the  molecule  parallel  and  perpendicular 
to  the  electric  field  of  the  laser  beam. 

Once  the  self  focusing  action  starts,  the  laser  beam  will 
have  to  travel  a certain  distance  through  the  medium  before  it 
will  come  to  focus.  It  can  be  shown  that  this  distance  is 
given  by  (22] 

P. 


->5 


(2.53) 


where 


and  P is  the  lower  of  the  optical  beam, 
o 


Let 


P = 2P  , , then 
o th 
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D = /2  D . 

o 

Now  if 

7 

= 10  V/cm 

n2  = 10-11  esu 

and  if  s = 1 cm  for  a CC>2  laser,  then 

D ■ 10  cm 
o 

so  that 

D = 14.14  cm. 

The  above  analysis  shows  that  self  focusing  would  be  very 
unimportant  for  semiconductor  materials,  where  the  thicknesses 
involved  are  only  a few  microns. 

It  should  be  noted  that  a beam  can  never  be  focused  to  a 
geometric  point  due  to  diffraction  effects. 

5.  PLASMA  INSTABILITIES  [27,  28] 

In  high  mobility  semiconductors  electrons  and  holes  can  be 
treated  as  free  particles  of  an  effective  mass  m* , interacting 
with  each  other  through  coulomb  forces.  Such  a system  can  be 
treated  like  a plasma  if  the  potential  energy  V due  to  the 
coulomb  interaction  and  the  kinetic  energy  K obey 

Zr  < 1 (2.54) 

Such  a condition  is  easily  satisfied  in  most  semiconductors 
and,  therefore,  semiconductors  can  be  assumed  to  be  containers  of 
electron  and  hole  plasmas.  A plasma  has  its  characteristic  fre- 
quency given  by 
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(2.5',) 


Whenever  the  frequency  of  an  interacting  electromagnetic  wave  is 
close  to  this  plasma  frequency,  several  different  types  of  in- 
stabilities can  be  excited  due  to  the  resonance  effect.  A brief 
description  of  two  of  the  important  instabilities  relevant  to  the 
present  work  is  given  below. 


Two  Stream  Instability  [29,  30] 

This  instability  has  been  treated  by  Nishikawa  [29]  and  also 
by  Kaw  and  Dawson  [30].  This  type  of  instability  can  occur  in  a 
plasma  in  which  the  ions  are  stationery  and  the  electrons  are 
moving  relative  to  it. 

The  basic  equations  of  motion  of  ions  and  electrons  are  the 
continuity  equation  and  the  momentum  conservation  relations  as 

iiven  below. 


— -■  + V • (n  . v . ) =0 
3t  3 3 


(2.5.  ) 


/ Dv  . \ 

n I — -1  + v . • Vv  . ) = - Vp  . + n .e  . (E  Sinujt  + E1 ) 
) 3 \ Jt  3 3 / 3 3 3 


- m .n  .v  .v  . 
3 3 3 3 


(2.57) 


where  Vp^.  gives  the  pressure  gradient,  E1  the  self  consistent 


field  of  the  electrons  and  the  ions  and  all  other  terms  have 
their  usual  meaning. 

Finally  E'  is  given  by  the  Poisson's  equation 


V • E'  = 4 tt V e ,n  . 

i 3 3 


(2.5H) 


Assuming  that  n^  (the  ion  density)  is  a slowly  varying  funet im, 
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of  time  and  neglecting 


one  can  solve  [30]  the  above  equations  to  get  the  following  fourth 
order  dispersion  relation  for  10 


1 = 


j;' 

k 

it 


m 


2 . ,2 
m.ui  U)  - kv  ) 
l o 


(2.59) 


where  v is  the  relative  velocity  of  electrons  w.r.t.  ions  and  k is 
o 

the  wave  vector  of  the  electromagnetic  wave.  If  all  the  roots 
of  Equation  2.59  are  real,  then  each  root  would  indicate  a pos- 
sible oscillation  given  by 
i (k  • r - uj  t) 

E’  oc  e n (2.60) 

If  some  of  the  roots  are  complex,  then  they  would  occur  in 

conjugate  pairs.  Let  these  roots  be  written  as 

U =6  + iy  (2.61) 

n n n 

The  time  dependence  of  the  E-field  will  now  be  given  as 
i (k  • r - a)  t)  y t 

E‘  <*  e n • e n (2.62) 

Equation  2.62  suggests  that  the  electrostatic  wave  can  grow 
exponentially  if  y^  is  positive.  This  would  give  rise  to  the  two 
stream  instability. 

Rewriting  Equation  2.59,  by  putting 

(A) 

x = — 

(A) 

P 
kv 

<i) 
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and 


y = 


o 


Figure  2-7.  Function  F(x,y)  in  the  two  stream  instability 

For  y = y0  > y,  the  line  F(x,y)  = 1 cuts  the  curve  at  four 
points  corresponding  to  four  real  roots.  However,  for  y = y^,  only 

't 

two  roots  are  real  and  the  other  two  must  be  complex,  giving  rise 
to  an  unstable  wave  in  the  plasma.  This  would  happen  for  small 
y,  i.e.  for  small  kvQ.  With  the  help  of  the  dispersion  relatipn 
it  can  be  shown  that  the  maximum  growth  rate  of  this  instability 
is  given  by 

y 

Oi 

> 

ft 
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Kaw  and  Dawson  [3]]  further  show  that  the  roots  of  the  dis- 
persion relation  are  complex,  i.e.  the  instability  occurs  if 


(2.65) 


where  T\  and  Tg  are  the  ion  and  electron  temperatures  and  can  be 
assumed  to  be  equal  in  case  of  semiconductors.  Hence  Equation 
2.65  can  be  written  as 


Convertinq  Equation  2.68  into 


power  threshold 
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th 


2 2 
m ui" 
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(2.66) 


(2.67) 
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(2.68) 


An  application  of  this  equation  to  calculate  the  damage 
threshold  for  CO^  lasers  and  for  Nd:YAG  laser  gives  results  that 
arc  2 to  3 orders  higher  than  the  observed  damage  thresholds. 


Parametric  Instability 

This  type  of  instability  was  originally  formulated  by  DuBois 
and  Goldman  [31,  32,  33,  34]  and  later  has  been  treated  by  several 
other  authors  [35,  36,  37], 

This  instability  arises  if  one  feeds  energy  to  the  plasma  at 

a rate  faster  than  it  can  dissipate  it.  Parametric  instability 

arises  when  an  intense  laser  beam  of  frequency  interacts  with 

the  plasma  of  resonance  frequency  cj  . If  in  - to  is  close  to  the 

R o R 

ion  acoustic  frequency  or  ion  optical  frequency  ojt  , the  ion  waves 

Lj 

will  absorb  energy  resonantly.  Similarly  ui  and  oj  can  interact 

Li 

to  pump  the  electron  plasma. 

Therefore,  monochromatic  pump  radiation  of  the  form 

! E [exp  i (k  • r - to  t)  + c • c] 

2 o o o 

will  modulate  the  plasma  at  the  beat  frequencies  to  ± [31]. 

The  longitudinal  polarization  in  a medium  under  the  action  of  the 
laser  field  will  be  a function  of  the  fields  at  to  + to  frequencies. 
The  modulation  of  the  longitudinal  polarization  can  be  neglected 
at  all  other  frequencies.  Therefore,  inside  the  medium 


-EL(k,J  = Pj  (k  ,„.) 


E.  (k  , to)  Xt  (k  »<o)  +E  • x^L  ( k 
L L o I,  o 


, (i)  ; k - k 

o o 


o 
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(2.70) 
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Similarly  one  writes 


P ( k,<j-u)  ) = -E  (k,o)-o)  ) = Et  (k,u)-w  )xT  (k,u)-u)  ) 

L o L oL  oL  o 

NL  -+ 

+ E *xT  (-k  ,-u)  ,k,u))E  (k,u))  (2.71) 

O L O O Lj 

Rewriting  Equations  2.70  and  2.71,  we  get 

->  ml 

e,  (k,w)ET  (k,w)+E  *XT  (k  ,w  ;k,u>-u>  )ET(k,w-w  ) = 0 (2.72) 

L L O L O O oL  O 

E •x!)L(-k  ,-w  ,k,w)E  (k,iii)+e  (k,w-w  )E  (k,w-w  ) = 0 (2.73) 

OLOO  L L OL  O 

where  E is  the  effective  longitudinal  electric  field  in  the  plas- 
L 

ma,  e (k,ui)  is  the  linear  longitudinal  dielectric  constant  and 
L 

NL 

XT  is  the  total  susceptibility  of  the  electron  ion  plasma.  xT 
L L 

is  the  nonlinear  longitudinal  susceptibility  that  connects  the 

longitudinal  polarization  to  E (u>  ± w ) . 

L O 

In  writing  Equations  2.72  and  2.73,  it  has  been  assumed  that 
— ► ->■  — ► 

inside  the  medium  k - k : k [31]. 

o 

Now  Equations  2.72  and  2.73  will  have  a non-trivial  solution 
if 


-*  ->  ^ml  > ->  -*NL  — 1 

i.  (k,w)-E  "X.  (k  ,u ) ;k,w-a)  )E  -xT  (-k  ,-io  ;k  ,w)  [eT  (k  ,u-oj  )]  =0 

L oLoo  ooLoo  L o 

(2.74) 

The  left  hand  side  of  Equation  2.74  is  also  called  the  non- 
linear longitudinal  dielectric  constant 

NL  1 
eL  (k  ,io)  . 

It  can  be  shown  that  [32] 
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then  the  term  in  the  square  brackets  in  Equation  2.75  can  be  shown 
to  be  equal  to 
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Substituting  this  in  Equation  2.76, 


e2 (k*E  )2  / k2\ 

EbL(k'“)  = L(k,<i,)  ' ' 2~f~  l“j)  tcL^,w-wo)] 

4m  w \ k ' 

e o 


e^L(k,w)  = e (k,w)  - — 
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(2.77) 


w \4  - I 

Cos  0 — —— 
wq  ) nckBT 


and  I is  the  incident  intensity, 
o 


Expression  2.77  will  be  used  in  Chapter  4 to  explain  C00 


damage  thresholds. 
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CHAPTER  3 


HALL  EFFECT  AND  SCATTERING  MECHANISMS 

Since  its  discovery  in  1879  [1] , Hall  effect  has  been  one  of 
the  most  powerful  methods  for  determining  the  carrier  mobility  and 
carrier  concentration  in  semiconductors.  Till  recently,  this  meth- 
od had  the  disadvantage  of  requiring  devices  of  specific  shape. 

But  in  1958,  Van  der  Pauw  [2]  extended  this  method  to  be  applied 
to  a sample  of  any  shape.  In  the  following,  we  shall  give  the 
theory  of  Hall  effect  and  Van  der  Pauw's  method  and  describe  the 
experimental  details  of  determining  mobility  and  carrier  concen- 
tration, before  and  after  laser  damage. 

1.  HALL  EFFECT  [3] 

When  a conductor  or  a semiconductor  carrying  current  is 
placed  in  a magnetic  field,  an  electromotive  force  is  produced 
across  it  in  a direction  perpendicular  to  both  the  current  and 
the  magnetic  field.  This  is  called  Hall  effect. 

Consider  an  n-type  semiconductor  bar  of  infinite  length  and 
width  d,  as  shown  in  Figure  3-1.  Let  be  the  current  density 
in  the  x direction  and  a magnetic  field  H be  applied  in  the  z- 
direction.  Let  n be  the  carrier  density. 

The  Lorentz  force  on  an  electron  in  the  bar  is  given  by 

F = e [E  + (v  x h) ] (3.1) 

where  E is  the  electric  field  and  v is  the  velocity  of  the  electron. 

As  no  current  can  flow  in  the  y-direction,  the  net  Lorentz 


force  in  this  direction  must  be  zero, 


Figure  3-1.  A current  carrying  semiconductor  bar 
of  infinite  length  placed  in  a 
magnetic  field 
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Equation  3.3  gives  the  magnitude  of  the  Hall  field.  Now  an 


electronic  flow  in  the  x-direction  sets  up  an  electric  field, 
given  by 

D = (3.4) 

where  o is  the  electrical  conductivity  of  the  semiconductor, 
or 

3^  = nep  E (3.5) 

ex 

From  Equations  3.3  and  3.5, 

= tan0  = -Hp  (3.6) 

e 

where  9 is  called  the  Hall  angle. 

Equation  3.3  is  usually  written  as 


where 

R = - — , and  is  called  Hall  constant.  In  deriving  the 
ne 

above  equation,  one  has  to  assume  that  the  energy  surfaces  of  the 
semiconductor  are  spherical,  which  is  rarely  true.  Therefore,  the 
expression  for  the  Hall  constant  should  be  corrected  as 


R = - — (3.8) 

ne 

where  r now  would  depend  upon  the  type  of  energy  surface  involved. 
In  general,  it  can  be  shown  [3]  that 


(3.9) 


where  t represents  the  statistical  average  of  the  electron  life- 


time ovor  the  proper  distribution  function. 


If  t^E  , which  is  true  in  acoustic  phonon  lattice  scatter- 
ing (E  being  the  energy) , then 

(t"  ) 2 3 

e 3u 


so  that 

R = (3. 

ne  8 

If  t as  in  case  of  ionized  impurity  scattering,  then 


(t  ) 2 
e 31o 


• R _ _ _i_  (315ii\ 

ne  \ 512  / 


From  Equations  3.10  and  3.11  it  can  be  seen  that  in  practice 
r is  very  close  to  unity  and  can  be  considered  so  in  experimental 
measurements . 

The  conductivity  of  an  n-type  sample  may  be  written  as 
o = -new  ( 3 . 1 : 


where  u is  called  conductivity  mobility, 
c 


Substituting  for  ne  in  Equation  3.8,  one  gets 


R = r • — 
a 

a UH 

or  p = R — = — 
c r r 


where  = R o and  is  called  Hall  mobility.  When  r - 1, 


|j  =-  p 
C H 


If  in  a sample  both  the  holes  and  electrons  are  present,  then 


it  can  be  shown  [5]  that 


1 Ptlh  ~ n>le 
e . , 2 

+ nMeJ 


(3.14) 


A determination  of  R experimentally  determines  the  carrier 
concentration  in  a sample.  If  the  conductivity  is  known,  then  one 
can  also  determine  the  mobility  of  the  carriers.  As  stated  in  the 
beginning,  this  method  requires  a sample  which  is  effectively  very 
long  as  compared  to  its  width.  Van  der  Pauw  has  extended  Hall  ef- 
fect method  in  such  a way  that  it  can  be  applied  to  a sample  of  an 
arbitrary  shape. 


2.  VAN  DER  PAUW  METHOD  [2] 

Consider  a sample  of  an  arbitrary  shape  as  shown  in  Figure 
3-2.  Let  d be  its  thickness  and  let  A,  B,  C and  D be  four  ohmic 
contacts  on  the  periphery  of  the  sample. 


Figure  3-2.  A sample  of  arbitrary  shape 
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Assume  the  following: 


(a)  The  contacts  are  at  the  circumference  of  the  sample. 

(b)  The  contacts  are  sufficiently  small. 

(c)  The  sample  is  homogeneous  in  thickness. 

(d)  The  surface  of  the  sample  is  singly  connected  and 
has  no  isolated  regions. 

Define 


AB , CD 


VC  VD 


then  one  can  show  by  using  elementary  electrostatics  that 


exp 


d 

d 

71  AB , CD  p 

+ exp 

^C.DA  p 

= 1 


(3.15) 


Here  p is  the  specific  resistance  of  the  sample.  Solving  Equa- 
tion 3.15  for  p, 


Trd  (RAB,CD  + ^C^a' 


5,  n2 


AB,CD 


*BC, 


DA 


where  f is  a function  of  the  ratio 


(3.16) 


AB  , CD 
, DA 

R 


only  and  satisfies 


AB , CD  V , DA  ,|  exp(f.n  2/f ) | 

= f arc  cosh  — 


RAB,CD  + ^C.DA  * 


If  RAB  CD  and  PBC  DA  are  nearly  e<3ual  (which  will  be  true  for 
a sample  of  nearly  square  shape) , f can  be  approximated  by 


1 - 


I(g.n2)  2 Un2)  3 I 
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If  RAB,CD  and  V,DA  are  eXaCtly  e9ua1'  then 
f = 1, 

and  from  Equation  3.16 

» - Vcd  <3a7> 

Now  if  a magnetic  field  is  applied  perpendicular  to  the 
sample,  then  the  equations 
div  ^ = 0 
0 

where  j is  the  current  density,  remain  valid.  As  the  contacts  are 
sufficiently  small  and  at  the  circumference  of  the  sample,  the 
outer  lines  of  flow  which  must  follow  the  circumference  fully  de- 
termine our  boundary  conditions.  Hence  the  lines  of  flow  do  not 
change  when  a magnetic  field  is  applied.  However,  the  magnetic 
field  causes  an  extra  potential  difference  between  any  two  arbi- 
trary points.  This  can  be  shown  from  Hall  effect  to  be 


curl  j 


AV  = 


VIIp 


g 

$ 

fe 


where  p is  the  Hall  mobility  and  I is  the  current. 

H 


UH  PH  I 


d_  AV 
> H 

d 


PH  " pH  ARBC,AC  (3.18) 

where  H is  the  magnetic  field  and  is  the  change  of  %d,AC 

due  to  the  magnetic  field. 

Using  Equations  3.17  and  3.18,  one  can  find  the  resistivity 
and  Hall  mobility  of  a sample,  provided  one  has 

(a)  A square  sample  of  the  semiconductor 

(b)  Four  ohmic  contacts  on  the  four  corners  of  the  sample. 


■ “ . ' ~ - ' w: 
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3.  FABRICATION  OF  HALL  DEVICE 


Bulk  semiconductor  materials  Si,  Ge,  GaAs  and  GaP  were  ob- 
tained commercially  in  wafers  ranging  from  0.01  cm  to  0.025  cm 

in  thickness.  GaAs  P,  was  obtained  in  epitaxial  (-300  pm)  wafer 
x 1-x 

form  on  intrinsic  (10®  S2cm)  GaAs  substrate. 

These  wafers  were  lapped,  polished  and  chemically  cleaned, 
using  procedures  developed  at  The  University  of  New  Mexico.  A 
description  of  these  procedures  is  given  elsewhere  [6] . 

Ohmic  contacts  were  evaporated  on  these  wafers  through  spe- 
cially prepared  masks  in  a vacuum  evaporation  system  at  pressures 
of  10”®  torr.  These  contacts  were  heat  treated  in  an  inert  atmos- 
phere from  5 to  15  minutes.  The  current-voltage  characteristics 
between  any  two  contacts  showed  good  ohmic  contacts. 

A diamond  subscriber  was  used  to  cut  the  wafers  into  (.25  x 
2 

.25)  cm  Hall  devices.  The  Hall  device  was  mounted  on  heat  con- 
ducting alumina  plates  by  using  heat  conducting,  electrically  in- 
sulating epoxy.  Four  strips  of  aluminum  were  evaporated  on  the 
four  corners  of  the  alumina  substrate.  Wire  bonds  were  made  be- 
tween these  four  plates  and  the  four  ohmic  contacts  on  the  sample. 
A typical  Hall  effect  device  is  shown  in  Figure  3-3. 


4. 


DESCRIPTION  OF  MOBILITY  APPARATUS 


The  Hall  effect  device  shown  in  Figure  3-3  was  mounted  on  a 
stage  specially  prepared  from  high  quality  oxygen  free  copper. 

This  stage  was  placed  in  perfect  thermal  contact  at  the  cold  tip 
of  the  displacer  of  a cryogenic  refrigerator,  as  shown  in  Figure 
3-4. 

The  displacer  is  surrounded  by  a vacuum  shroud  that  is  evac- 
-4 

uated  to  10  torr  by  using  a mechanical-diffusion  pump  system. 

Four  wires  are  fed  through  to  make  electrical  connections  to  the 
Hall  device.  A heater  element  is  wound  over  the  cold  tip  to  con- 
trol its  temperature,  which  can  be  adjusted  from  12 °K  to  300 °K. 

A chromel  vs.  gold — 0.07  atomic  percent  iron  thermocouple  is  used 
to  measure  the  temperature  of  the  cold  end  to  an  accuracy  of  ±0.2%. 

The  displacer  assembly  is  mounted  between  adjustable  pole 
pieces  of  a water  cooled  electromagnet.  The  magnet  is  capable  of 
giving  fields  up  to  2.5  teslas. 

Resistivity  and  mobility  were  determined  using  the  following 
steps . 

Step  1:  A current  of  a few  milliamperes  was  passed  through 

two  adjacent  contacts  of  the  sample  by  using  dc  batteries.  The 
Hall  induced  electromotive  force  was  measured  across  the  other 
two  contacts.  The  current  was  then  reversed  in  direction  and  Hall 
voltage  was  remeasured.  This  was  repeated  for  all  combinations 
of  adjacent  contacts  and  the  average  of  the  eight  readings  was 
taken  to  give  of  Equation  3.17  and  hence  to  determine 

resisti vi ty  p . 

step  2:  The  current  is  now  passed  through  two  diagonal  con- 


48 


Device 


tacts  of  the  sample  and  the  Hall  voltage  is  measured  across  the 
other  two  diagonal  contacts.  Then  a magnetic  field  is  applied  in 
a direction  perpendicular  to  the  direction  of  current  and  the 
change  in  Hall  voltage  determined.  The  magnetic  field  is  reversed 
and  the  change  in  Hall  voltage  remeasured. 

Step  3:  The  current  through  the  device  is  reversed  and  Step 

2 is  repeated. 

Step  4:  The  current  is  now  passed  through  the  other  two 

diagonal  contacts  and  Steps  2 and  3 are  followed. 

Average  of  these  eight  readings  is  taken  to  give  AR^  A(,  of 
Equation  3.18  and  thus  Hall  mobility  is  determined. 

5.  SCATTERING  MECHANISMS 

The  mobility  of  carriers  in  a semiconductor  is  determined 
by  the  scattering  mechanisms  in  the  lattice.  Radiation  damages 
the  lattice  structure  of  semiconductors  due  to  which  electronic 
properties  are  degraded.  Particle  damage  is  known  to  cause  change 
in  carrier  concentration,  lifetime  and  mobility  of  a semiconductor 
due  to  the  creation  of  interstitial-vacancy  pairs  or  complexes. 
These  changes  may  be  significant  depending  on  the  amount  of  ra- 
diation dose.  The  defects  introduced  may  or  may  not  dominate 
the  electrical  properties  at  room  temperature  because  of  the  lat- 
tice vibrations.  However,  at  low  temperatures  lattice  vibrations 
can  be  killed  and  the  effect  of  the  introduced  defects  easily 
observed.  Therefore,  in  this  chapter  the  temperature  dependence 
of  different  scattering  mechanisms  that  control  the  carrier  mo- 
bility and  electrical  conductivity  of  the  semiconductor  lattice 


will  be  examined.  This  dependence  is  what  eventually  determines 
the  nature  of  radiation  induced  defects. 

Scattering  mechanisms  can  be  divided  into  three  parts. 

(a)  Lattice  scattering  [7,  8,  9] 

(b)  Impurity  scattering:  ionized  [10,  11]  and  un-ionized  [12] 

(c)  Defects  [13,  14] , intervalley  [15]  and  carrier-carrier 
scattering  [16] 

( a)  Lattice  Scattering 

The  vibration  spectrum  of  semiconductors  has  two  branches, 
namely  optical  and  acoustical.  At  moderate  and  low  temperature, 
the  thermal  energy  available  is  insufficient  to  excite  the  high  fre- 
quency optical  mode  and,  therefore,  the  scattering  of  carriers  can 
be  due  to  acoustical  vibrations  only.  One  considers  only  the  scat- 
tering of  charged  carriers  by  longitudinal  mode  vibrations,  because 
the  effect  of  transverse  modes  can  be  shown  to  be  negligible. 

The  passage  of  a longitudinal  vibration  through  a crystal 
would  give  rise  to  alternate  regions  of  compression  and  extension. 

A compressed  region  would  exhibit  an  increase  in  band  gap,  while 
an  extension  would  cause  a decrease  in  the  band  gap.  This  is 
shown  in  Figure  3-5. 

Along  the  length  of  the  crystal,  therefore,  the  band  gap  as 
a function  of  distance  would  look  like  that  shown  in  Figure  3-6 (a) 
and  can  be  approximated  by  Figure  3-6  (b)  . One  may  then  calculate 
the  reflection  probability  for  an  electron  incident  upon  a single 
step  of  height  6Ec  as  shown  in  Figure  3-6 (b).  Neglecting  Doppler 
effect  due  to  a barrier  moving  at  the  speed  of  sound,  it  can  be 
shown  from  elementary  quantum  mechanics  that 
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ft 


E being  the  energy  of  the  incident  electron.  If  the  step  6E(  is 
small,  which  can  safely  be  assumed,  then  k^  = kQ.  With  this  ap- 
proximation, Eguation  3.19  becomes 


Now  to  a first  approximation 


(3.22) 


(3.23) 


where 

6V 

V 

o 

represents  compressional  or  extensional  strain  and  n*  is  the  de- 
formation potential  constant.  If  6p  is  the  maximum  pressure 
created  by  the  compression  or  extension,  then  the  stored  strain 
energy  is 


<5E  = \ <5p6v  = A'kT  (3.24) 

2.  B 

where  A*  is  a constant.  The  term  X'k^T  arises  because  the  source 
of  this  strain  energy  is  thermal.  Also  by  definition  the  bulk 
compressibility  is  given  as 


B = 


1 6v 

V 6p 
o 


From  Equations  3.24  and  3.25,  one  has 

( £)'  - 


>2  2 B A 'k  T 

B 


Substituting  this  in  Equation  3.23  gives 


(3.25) 


(3.26) 
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(3.27) 
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From  Equations  3.22  and  3.27  one  has 


V2n2BA'kBT 
R 2 — 

2V  (fi  k ) 
o o 


(3.28) 


If  i/2  is  the  linear  dimension  of  the  volume  V , Xg  the 
mean  free  path,  then  in  going  through  a distance  i/2,  the  scat- 
tering probability  is  i/2\^.  Therefore,  replacing  R in  Equation 
3.28  by  S./2X produces 

t me*2n,2BX'kBT 

2X~  _ j-4.  4. 

e V (h  k ) 
o o 


or  X = 
e 


8m*2X'6k„Tn’2 
e B 


(3.29) 


where  we  have  replaced  kQ  by  2i\/i  and  by  i /8  . 

The  collisions  between  electrons  and  phonons  are  very  simi- 
lar to  perfectly  elastic  collisions  between  two  hard  spheres  and 
can  be  described  by  a velocity  independent  mean  free  path.  As- 
suming the  velocity  distribution  of  electrons  is  Maxwell  Boltzmann 
it  can  be  easily  shown  that  the  mean  free  time  is  given  by  [18] 


(3.30) 


where  c is  the  mean  thermal  speed  of  the  electron.  From  Equations 
3.29  and  3.30  one  has 


6 ' 3/8^  m*3/?X-B(k  T)3/2n.2 
e B 


(3.31) 


Even  though  the  above  equation  has  been  derived  qualitatively, 
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it  is  not  very  much  different  from  the  quantum-mechanical  ex- 
pression [19],  which  is  given  as 


t = 
e 


/r 


fi4C 


£11 


m*3/2  (k0T) 3/2n'2 


(3. 32) 


where  is  the  elastic  constant  for  a longitudinal  extension  in 
the  [110]  direction.  The  corresponding  mobility  is 


et 


/r 


eR4C 


ll 


m* 

e 


5/2 


(kBT) 


3/2  ,2 


(3.33) 


- 3/2 

The  theoretical  prediction  ueaT  is  only  approximately 


true  and  in  actual  practice  the  variation  is  found  to  be  stronger 
in  Ge  and  Si.  This  can  be  attributed  to  the  following  reasons. 

(i)  In  p-type  Ge  and  Si,  optical  scattering  cannot  be 
ignored. 

(ii)  In  n-type  Ge  and  Si,  energy  surfaces  are  ellipsoidal 
and  effects  of  multivalley  scattering  cannot  be 
ignored  [20]. 


(b)  Impurity  Scattering 

Ionized  Impurities.  The  determination  of  the  effect  of 
ionized  impurity  scattering  on  mobility  of  electrons  is  based 
on  the  theory  of  scattering  of  charged  particles  by  the  coulomb 
potential  of  the  nucleus,  which  was  originally  developed  by 
Rutherford  [21]  . He  showed  that  the  scattering  cross  section  of 
a charged  particle  of  charge  e by  a fixed  nucleus  of  charge  Ze 
is  given  by 
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Figure  3-7.  Rutherford  scattering  of  an  electron 
from  a nucleus  of  charge  Ze 

The  number  of  collisions  an  electron  can  make  per  unit  time 

into  a solid  angle  d&  can  be  shown  to  be  Nv  o(0)dft  where  N is  the 

o 

number  of  charged  scattering  centers  per  unit  volume.  In  each 
such  collision  the  momentum  change  is  only  mvQ(l  “ Cos0),  as  can 
be  seen  from  Figure  3-7.  The  relaxation  time,  however,  is  the 
time  during  which  the  average  forward  momentum  is  zero.  Hence  in 
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this  sense,  the  effective  number  of  collisions  is  given  by 

d I 


dn  ,,  = Nv  o(0)  (1  - Cos0)d£2 
ef  f o 


'(t) 


(3.36) 


Equation  3.36  gives  the  differential  of  inverse  time  period  for 
electrons  of  initial  velocity  vq  scattered  through  an  angle  0 into 
a solid  angle  dft.  The  total  relaxation  time  may  be  found  by  inte- 
grating Equation  3.36  over  the  solid  angle  and  averaging  over 
velocities . 

Assuming  the  scattering  to  be  independent  of  the  azimuthal 
angle,  integration  of  Equation  3.36  over  the  solid  angle  from  polar 
angle  0 to  u gives 


, * 

= 2ttNv  / o(0)  (1  - Cos0)Sin0dO 

T ( V ) O 1 


Jf. 
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fe 


= 8ttNv 
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o|  , „2 
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em*v 


1 + 


2Ze2N  1/3 ^ 


(3.37) 


The  lower  limit  of  integration  has  not  been  taken  to  be  zero, 
because  in  a crystal  the  coulomb  potential  of  an  ion  has  to  be  cut 

off  at  a certain  distance  from  it  due  to  the  presence  of  other  ions. 

-1/3 

It  is  reasonable  to  take  this  distance  to  be  d/2  where  d = N 
is  the  mean  distance  between  impurity  atoms. 
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and  accordingly. 
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This  expression  has  already  been  used  to  get  Equation  3.37. 

It  is  now  necessary  to  average  Equation  3.37  over  the  Max- 
wellian distribution  of  velocities. 
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The  integral  in  Equation  3.38  can  be  solved  only  approximately. 


This  is  done  by  assuming  that  the  logarithmic  term  varies  much  more 


slowly  with  vq  as  compared  to  the  vq  term.  Making  this  approxi- 


mation 
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(3.39) 


This  formula  is  called  Conwell-Weiskopf  formula  [10] . Equa- 


tion 3.39  shows  that  the  mobility  due  to  ionized  impurity  scatter- 

3/2 


ing  is  directly  proportional  to  T and  inversely  proportional  to 
impurity  concentration.  Therefore,  for  low  impurity  concentrations, 
only  lattice  scattering  will  be  dominant. 
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Neutral  Impurities.  The  scattering  by  neutral  impurities  is 
analogous  to  the  scattering  of  electrons  by  hydrogen  atoms. 
Erginsoy  [12]  has  shown  that  the  relaxation  time  for  this  process 
is  given  by 


-L  = / 201jlL 
t I 2 2 In ' 

e \ e m*  I 

e ' 


where  N is  the  number  of  neutral  impurities, 
n 


Hence  mobility  dependence  due  to  un-ionized  impurities  is 
given  by 


V*  = 


ei  3m* 

e e e 


20R  eN 


(3.40) 


This  type  of  scattering  is  important  only  at  very  low  tem- 
peratures where  most  of  the  impurities  are  un-ionized. 


6 


( c ) Defect  Scattering,  Intervalley  Scattering  and  Carrier- 

Carrier  Scattering 

Defects  such  as  dislocations  scatter  carriers  as  a result  of 
the  strain  around  the  dislocation  which  gives  rise  to  a deforma- 
tion potential.  This  has  been  treated  by  Dexter  and  Seitz  [13]. 
The  predicted  temperature  dependence  of  the  mobility  arising  from 
such  dislocations  is  linear  in  T. 

Dislocations  also  act  as  acceptors  and,  therefore,  act  as 
negatively  charged  centers  causing  additional  scattering. 

Intervalley  scattering  involves  the  inelastic  transfer  of 
an  electron  from  one  band  at  a certain  crystal  momentum  to  the 
part  of  the  same  band  at  a different  crystal  momentum.  This  has 
been  further  discussed  by  C.  Herring  [15] . 

Carrier-carrier  scattering  can  be  neglected  under  ordinary 
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circumstances.  However,  carrier-carrier  collisions  randomize  the 
electron  velocity  distribution,  thereby  changing  the  probability 
of  electron  scattering  through  other  mechanisms.  Such  a treatment 
has  been  given  by  Spitzer  and  Harm  [16]. 
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CHAPTER  4 


RESULTS:  CW  CC>2  LASER  DAMAGE 

The  experimental  results  on  damage  thresholds  and  the  changes 
in  electrical  properties  of  the  semiconductors  due  to  C02  laser 
radiation  are  given  in  this  chapter.  The  results  obtained  are 
explained  theoretically. 

The  damage  thresholds  were  obtained  by  radiating  the  chemi- 
cally cleaned  semiconductor  samples  in  the  form  of  Hall  devices. 
The  damage  criteria  used  have  been  described  in  Chapter  2.  The 
CW  CO  laser  beam  used  was  Gaussian  in  space  with  a diameter  of 
12  mm.  A salt  lens  of  25  cm  focal  length  was  used  to  concentrate 
the  beam  on  the  sample.  The  sample  was  displaced  from  the  focus 
so  that  the  whole  area  could  be  uniformly  radiated.  The  experi- 
mental arrangement  is  shown  in  Figure  4-1. 

The  absorption  coefficients  [1]  for  most  semiconductors  at 
10.6  iim  art;  very  small  and,  therefore,  it  can  be  assumed  that  the 
energy  is  absorbed  by  the  bulk  uniformly. 

1.  DAMAGE  THRESHOLDS 

The  experimentally  observed  damage  thresholds  for  different 
semiconductors  are  listed  in  Table  4-1.  Also  given  in  the  same 
table  are  their  pre-damage  mobility,  carrier  concentration  and 
band  gap. 

From  this  table  we  observe  that  the  damage  thresholds  are 
ill  of  the  same  order  and  are  quite  insensitive  to  the  original 
arrier  concentration  and  mobility. 
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Figure  4-1.  Experimental  arrangement  used  to  radiate  the  sample 


ELECTRICAL  PROPERTIES 


Conductivity 

The  electrical  conductivity  of  Ge  as  a function  of  temperature 
before  and  after  laser  damage  is  shown  in  Figure  4-2.  It  is  seen 
that  there  is  a significant  decrease  in  conductivity  due  to  laser 
damage , particularly  at  low  temperatures.  This  shows  that  the  de- 
fects introduced  by  damage  are  very  effective  scattering  centers  at 
lower  temperatures.  The  lattice  vibrations  dominate  the  created 
defects  at  room  temperature. 

Carrier  Concentration 

The  carrier  concentration  versus  temperature  before  and  after 
laser  damage  is  plotted  for  Si  and  Ge  in  Figures  4-3  and  4-4  respec- 
tively. A magnified  view  of  Figure  4-3  is  given  in  Figure  4-5  to 
show  the  significant  carrier  removal  due  to  laser  damage. 

In  general,  the  carrier  concentration  is  decreased  both  in 
Ge  and  Si  at  all  temperatures.  This  shows  that  the  laser  created 
scattering  centers  give  rise  to  localized  energy  levels  in  the  for- 
bidden band  which  act  as  carrier  traps. 

Mobility 

The  mobility  changes  due  to  laser  damage  were  very  small  at 
room  temperature.  However,  as  the  temperature  is  lowered,  a sig- 
nificant degradation  of  mobility  occurs  both  in  Ge  and  Si  (see 
1 igures  4-6  and  4-7) . 

This  confirms  the  idea  that  the  introduced  defects  are  more 
effective  scattering  centers  at  lower  temperatures. 
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Figure  4-3.  Carrier  concentration  versus  temperature  for  Si  be- 
fore and  after  laser  damage  (The  sample  was  radiated 
at  0.96  x 10^  watts/cm2.  The  original  carrier  con- 
centration was  ~10l5  at  room  temperature.) 
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Figure  4- 


Electron  mobility  versus  temperature  for  Si  before 
and  after  damage  (Other  parameters  are  same  as  in 
Figure  4-3.) 
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3. 


THEORETICAL  EXPLANATION  OF  DAMAGE  THRESHOLDS  [2] 


i 

i 


i 

t 


A systematic  analysis  of  different  damage  mechanisms  given  in 
Chapter  2 shows  the  following. 

(a)  All  band  to  band  transitions  are  energetically  forbidden, 
because  the  photon  energy  for  CO  2 laser  is  only  0.117  eV, 
which  is  much  smaller  than  any  of  the  band  gaps  involved. 
Transitions  to  surface  states  are  possible,  however. 

(b)  Simple  Joule  heating  is  insufficient  to  raise  the  tem- 
perature of  the  sample  to  melting  point  in  a time  of  about 
0.1  second  for  which  the  samples  were  irradiated. 

(c)  Multiphoton  absorption  cross  sections  are  too  low  [3] 
to  cause  any  significant  absorption. 

(d)  Damage  to  doped  and  undoped  samples  occurs  at  the  same 
power  levels.  This  rules  out  the  occurrence  of  damage 
due  to  free  carriers.  Free  carrier  absorption,  however, 
can  be  a cause  of  raising  the  electron  temperature  in  the 
early  stages  of  damage  process. 

(e)  Lucky  electron  theory  of  avalanche  cannot  apply  to  high 
mobility  semiconductors,  because  of  the  fact  that 

u)ie  >>  1.  This  inequality  rules  out  the  possibility  of 
a hard  momentum- reversing  collision  during  half  period 
of  the  laser,  essential  for  the  buildup  of  an  electron 
avalanche.  Shatas  et  al.  [4]  have  also  ruled  out  elec- 
tron avalanche  as  a probable  damage  mechanism  in  GaAs 
and  NaCl  crystals. 

This  leaves  one  with  plasma  instabilities.  Upon  examination 
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of  several  different  types  of  plasma  instabilities,  it  was  deter- 
mined that  the  theory  of  parametric  instability  of  DuBois  and  Gold- 
man [5] , when  combined  with  the  theory  of  anomalous  absorption  of 


Dawson  and  Oberman  [6] , explains  the  observed  damage  thresholds 
satisfactorily.  ( 

This  kind  of  instability  can  be  excited  in  a system  if  the 
energy  is  fed  into  the  solid-state  plasma  at  a rate  faster  than  the 
system  can  dissipate  it.  The  intense  laser  beam  of  frequency  <d 
(2.83  x 1013  Hz  in  this  case)  interacts  nonlinearly  with  the  elec- 
tron plasma,  that  has  a resonance  frequency  given  by 

"r  ’ (“e  * 3kV)'‘ 

(where  u>e  is  the  Langmuir  frequency  of  the  electron  plasma  and  se 

is  the  thermal  velocity  of  the  electron)  and  pumps  the  ion  plasma 

at  the  beat  frequency  (w  - m ) . If  w and  u>_  are  very  nearly 

oR  o R 

equal  , then  ( co  - oi  ) is  small  and  is  characteristic  of  the  ion 
o R 

acoustic  frequency.  Similarly,  the  laser  beam  can  interact  with 
the  ion  acoustic  wave  and  pump  the  electron  plasma.  If  this  type 
of  instability  has  a positive  growth  rate,  then  it  can  cause  anom- 
alous absorption  and  heating  of  the  sample  with  subsequent  lattice 
damage . 

The  calculation  of  the  damage  thresholds  is  based  on  finding 
the  nonlinear  susceptibility  in  terms  of  the  laser  beam  parameters 
(Eq,  kQ)  and  electron  plasma  parameters  (n , T,  we)  . 

It  was  shown  in  Chapter  2 (see  Equation  2.77) 
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for  k <<  kD,  A <1  where  eL(k,io)  is  the  linear  longitudinal  dielec- 

NL  -* 

trie  constant  of  the  electron  plasma  and  e (k,aj)  is  the  nonlinear 
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dielectric  constant  under  the  action  of  the  intense  laser  beam.  Here 
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in  which  is  the  usual  Langmuir  frequency  of  the  plasma  and  is 
the  frequency  of  the  incident  laser  field  and  Iq  its  density,  n 
is  the  electron  concentration  of  the  plasma.  Also,  kQ  is  the  in- 
verse of  the  Debye-Huckel  shielding  length. 

The  nonlinear  dielectric  constant  has  two  parts.  The  real 
part,  when  equated  to  zero,  would  give  the  resonance  frequency  of 

the  plasma  and  is  always  very  close  to  ojq.  The  imaginary  part 

, NL 

gives  the  damping  rate  (Y  /uj^) 
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(4.3) 


where 


Im  ( e (k,wT  ))  = — 

\ L LI  a) 

P 


(4.4) 


and  is  ordinary  Landau  damping.  In  the  case  of  a solid  state  plasma, 
it  can  be  approximated  as 
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P p I>C 

where  is  the  collision  time  of  the 
phonons . 


(4.5) 

carriers  with  the  lattice 


Im 
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can  be  estimated  by 


ft)  ft) 


where  Y^/ar  and  Y^/w  give  the  damping  rates  for  acoustic  and 
optical  phonons,  respectively  [7]. 

Hence 


This  damping  rate  is  negative  (i.e.,  we  have  the  growth  of 
instability)  if 


or  writing  in  terms  of  [xswer  threshold 


(4.6) 


Now  gj  .i  and  for  this  to  be  true  the  plasma  concentration 

op  1 

should  be  of  the  order  of  1018  to  1019  cm-3.  But  the  starting 
carrier  concentrations  here  are  much  less  than  this.  Therefore,  a 
new  excitation  mechanism  is  proposed  which  would  alleviate  enough 
electrons  from  the  valence  band  to  the  conduction  band  via  a con- 
tinuous distribution  of  surface  states  within  the  band  gap.  Tins 
will  be  disucssed  at  the  end  of  this  section. 

In  a solid  state  plasma  Drude  absorption  causes  the  electron 
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Table  4-2. 

Theoretical  and  Experimental 
for  Different  Semiconductors 

Damage  Thresholds 

Calculated 

Experimental 

Damage  Thresholds 

Damage  Thresholds 

Material 

-2 

in  W cm 

-2 

in  W cm 

Si  (n-type) 

4 

7.2  x 10 

1.1  x 104 

Si  (p-type) 

9.8  x 104 

2.2  x 104 

Ge 

4.08  x 104 

8.5  x 103 

GdAs 

7. 14  x 103 

6 x 103 

GaAsP 

5.7  x 103 

3.5  x 103 

Gap 

8 x 104 

1.33  x 103 

KC1 

4 5 

10  -►  10 

4 5 

10  + 10 

3 tLH  ’ 


Photon  Induced  Excitation 


It  is  well  known  that  the  width  of  the  band  gap  for  a semi- 
ductor  is  much  less  near  the  surface  as  compared  to  within  the  bulk, 
e.g.,  the  band  gap  for  GaAs  near  the  surface  is  0.30  eV  to  0.36  eV 
and  1.6  eV  within  the  bulk  [9J.  Also  within  the  band  gap,  there  is 
a continuous  distribution  of  surface  states  having  a density  of  about 
1012  to  1014  cm  2 (eV)  The  lifetimes  for  these  states  are  3 to 

4 orders  longer  as  compared  to  for  C02  laser.  The  electrons  in 
the  states  near  the  Fermi  level  in  an  n-type  semiconductor  are  al- 
most resonantly  excited  to  the  conduction  band.  The  emptied  states 
are  filled  again  by  the  electrons  excited  by  the  photons  from  the 
states  lying  below  the  Fermi  level.  This  streaming  process  can  pro- 
duce electron  concentrations  of  10'*'®  to  10*"®  cm  2 within  milli- 
seconds or  less.  The  quantum  yields  and  transition  probabilities 
for  such  processes  have  been  calculated  by  Kane  [10]  and  experi- 
mentally verified  by  Gobeli  and  Allen  [11].  The  quantum  yield  has 

been  shown  to  be  of  the  order  (hv-E  )n  where  n = 1 to  5/2,  E is 

s s 

the  separation  between  two  surface  states  and  hv  is  the  photon 
energy.  This  high  quantum  yield  makes  the  proposed  mechanism  highly 
probable . 

4.  EXPLANATION  OF  ELECTRICAL  DAMAGE  [12] 

Radiation  damage  in  semiconductors  produces  lattice  defects 
in  the  form  of  the  interstitials,  vacancies,  di vacancies  or  va- 
cancy-impurity complexes  depending  upon  the  nature  of  the  radia- 
tion and  the  semiconductor  involved.  In  case  of  neutron,  electron 
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and  y-ray  damage  of  semiconductors,  the  creation  of  these  defects 


can  easily  be  explained  by  considering  the  collision  between  the 
individual  radiation  particle  and  the  lattice  atoms.  In  the 
CW  CO 2 laser  damage,  however,  a single  photon  has  only  .117  eV  of 
energy  and,  therefore,  cannot  cause  displacement  of  atoms.  As 
discussed  in  the  previous  section,  the  damage  mechanism  in  the 
present  case  is  the  anomalous  absorption  pf  laser  energy  by  the 
solid-state  electron  plasma.  The  anomalously  heated  plasma  gives 
up  its  energy  to  the  lattice  atoms  primarily  through  Landau  damp- 
ing [13] . At  the  given  damage  thresholds  the  energy  gained  by 

' 19  3 

the  lattice  is  enough  to  displace  10  atoms/cm  , assuming  that 
the  displacement  energy  per  atom  is  25  eV  [14] . Most  of  the  dis- 
placed atoms  end  up  in  other  normal  lattice  sites,  while  about  1 
in  106  finally  occupies  an  interstitial  site,  thereby  creating 
interstitial-vacancy  pairs,  some  di vacancies  or  vacancy-impurity 
complexes  may  also  be  produced  in  the  process.  The  net  result  of 
the  formation  of  these  defects  is  the  creation  of  localized  elec- 
tronic levels  in  the  bandgap  that  act  as  trapping  centers. 

Earlier  ideas  concerning  the  energy  levels  associated  with 
the  interstitial  and  the  vacancy  were  based  on  their  single  ioniza- 
tion energies.  However,  in  1951,  it  was  shown  by  James  and  Lark- 

Horovitz  [15]  that  the  second  ionization  energies  of  the  inter- 
stitial and  vacancy  also  play  a dominant  role  in  determining  the 
activation  energies  of  the  localized  levels  associated  with  the 

»„ 

radiation  damage.  Based  on  these  arguments,  the  following  local- 

*■ 

ized  levels  can  be  predicted: 
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(a)  A shallow  vacant  level  corresponding  to  the  first 
ionization  energy  of  Si  or  Ge,  located  at  approxi- 
mately 0.05  eV  below  the  conduction  band. 

(b)  A deep,  vacant  level  corresponding  to  the  second 
ionization  energy  of  the  interstitial  or  vacancy 
in  Si  or  Ge . 

The  two  remaining  levels  lie  very  close  to  the  valence  band 
and  shall  always  be  occupied.  Predicting  the  precise  location 
of  these  levels  with  the  help  of  theoretical  calculations  is  ex- 
tremely difficult  because  of  the  complex  nature  of  the  wave  func- 
tions associated  with  the  defects  and  their  large  perturbing 
effect.  In  the  discussion  that  follows,  we  shall  assume  that 
only  the  empty  levels  can  cause  any  carrier  removal  through  the 
equation 

electron  + trap  t occupied  trap. 

Le  n be  the  electron  concentration  and  N be  the  concentra- 
tion of  the  deep  traps.  N is  also  the  concentration  of  the  shal- 
low traps.  At  a given  temperature,  only  a certain  number  of 
shallow  traps  are  occupied,  while  all  the  deep  traps  are  occu- 
pied at  all  temperatures.  Let  Nfc  be  the  concentration  of  the 
occupied  traps,  then 

n + (N  - Nfc)  ^ Nt 

electron  vacant  traps  occupied  traps  ’ (4.8) 

Applying  the  law  of  mass  action  to  this  reaction,  one  has 


Kt  n(n  - Nfc) 
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or 


Nt  1 + nK, 


where  Kfc  is  a constant  and  can  be  shown  to  be  [13] 


h exp[Es/kBT] 

(16Ttm*k  T)  3/2 
e b 


(4.11) 


where 


Eg  = Depth  of  the  shallow  traps  below  the  conduction  band 

T = Temperature  in  absolute  degrees 

k0  = Boltzmann  constant 

m*  = Electron  effective  mass, 
e 

Now  at  any  temperature , the  carrier  concentration  can  be 


written  as 


n = n - N.  - N (4 

o t 

where  n is  the  original  electron  concentration, 
o 

Substituting  for  Nfc  from  Equation  4.10  in  Equation  4.12 


n = (n  -2N)K  - 

2K^  o t 


i +V(  (nQ-2N) Kt  - 1)  + 4K 


:t(no-N) 


Equation  4.13  can  be  used  to  predict  carrier  concentration 
versus  temperature,  provided  N is  known.  N can  be  found  from 
the  mobility  data  as  follows. 

As  discussed  earlier,  the  Frankel  defects  act  like  charged 
scattering  centers.  Therefore,  one  can  apply  Brook-Herring  [16] 
theory  according  to  which  the  mobility  due  to  ionized  impurities 


» 


is  given  by 


UI  2N  B 


(4.14) 


27^  x (,)-a/;(V)3/2e2 

Z2e3m*>5 

e 


and  e = Dielectric  constant  of  the  semiconductor 


Z = Nuclear  charge  of  the  impurity 


(4.15) 


e = Electronic  charge 
B = In (1  + b)  - b/(l  + b) 


6 (kgT)  m* 
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(4.16) 


(4.17) 


At  any  temperature,  lattice  vibrations  also  contribute  to 
scattering  and  the  total  mobility  can  be  written  as 


(4.18) 


where  is  the  mobility  as  determined  by  lattice  vibrations. 

The  radiation  damage  does  not  affect  the  lattice  scattering, 
therefore,  any  change  in  total  mobility  can  be  brought  about  only 
by  changing  the  total  number  of  charged  scattering  centers. 


or  from  Equation  4.14, 


.*Z-  at&* 


r — 1 i 

ANj  gives  the  total  number  N of  Frankel  defects  created  due  to 
the  laser  damage 

N-A4(i)  (4.19) 

Using  Equation  4.19  and  the  mobility  versus  temperature 
curves  at  several  temperatures,  mean  values  of  N for  Ge  and  Si 
were  determined  to  be  3.32  x 10^  cm  ^ and  8.1  x 10^  cm  ^ , re- 
spectively. Using  these  values  of  N and  Equation  4.13,  it  is 
found  that  the  predicted  n versus  T dependence  coincides  with  the 
experimental  for  Eg  = 0.18  eV  in  the  case  of  Ge  and  for  Eg  = 0.36 
eV  in  the  case  of  Si  (see  Figures  4-3  and  4-4) . 

We  see  that  the  experimentally  determined  trap  depth  of 
0.18  eV  differs  considerably  from  the  theoretically  predicted 
trap  depth  of  0.05  eV.  However,  theoretical  prediction  is  based 
on  an  extremely  simple  hydrogenic  model  which  is  not  true  in  ac- 
tual practice.  Neutron  [17]  and  y-ray  [18]  damage  studies  in  Ge 
have  also  shown  the  existence  of  this  localized  level  at  about 
the  same  depth.  Similar  radiation  damage  studies  in  the  case  of 
n-type  Si  have  shown  [19,  20]  the  existence  of  two  traps:  one 
lying  at  0.16  eV  and  the  other  at  0.40  eV  below  the  conduction 
band,  the  former  being  produced  much  more  frequently  than  the 
latter.  Moreover,  D.  E.  Hill  [21]  has  found  evidence  for  a local- 
ized state  at  0.36  eV  below  the  conduction  band,  which  coincides 
with  the  present  study. 
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CHAPTER  5 


RESULTS:  Nd:YAG  LASER  DAMAGE 


This  chapter  presents  the  results  on  laser  damage  of  semi- 
conductors at  1.06  pm  wavelength.  The  beam  divergence  was 
-4 

2 x 10  radians  and  its  duration  between  half  power  points  was 

-35  ns.  The  pulse  shape  was  Gaussian  in  space  and  time,  as  shown 
in  Figure  5-1. 


Figure  5-1.  Nd : YAG  Laser  Pulse  (width  ~35  ns) 

The  manner  in  which  the  samples  were  irradiated  is  shown  in 
Figure  5-2.  Different  neutral  filters  were  used  in  the  path  of 
the  beam  to  adjust  the  irradiance  at  the  sample.  A reflected  beam 
off  a beam  splitter  was  fed  to  an  energy  meter  to  measure  the 
energy  incident  on  the  sample.  The  criteria  used  for  damage  are 
the  same  as  in  References  (1,  2]. 

A lens  of  focal  length  50  cm  was  used  to  focus  the  beam  on 
the  sample.  The  size  of  the  focal  spot  on  the  sample  is  calcu- 
lated from  [3] 

s = f • 0 (5.1) 
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Figure  5-2.  Experimental  arrangement  used 
to  radiate  the  samples 


where  f is  the  focal  length  of  the  lens  and  0,  the  beam  divergence 
Equation  5.1  is  true  only  at  very  small  distances  from  the  laser. 
This  condition  was  satisfied  in  the  experimental  arrangement.  Dif 
ferent  semiconductors  with  their  solid  state  parameters  and  their 
respective  experimental  damage  thresholds  are  shown  in  Table  5-1 . 


Table  5-1.  Damage  threshold  for  different 
semiconductors  at  1.06  microns 


Semiconductor 

Mobility 
(cm1 2 * * * * *(V/sec)  ■*■) 

Carrier 

concentration 

. -3, 

(cm  ) 

Band 

gap 

(eV) 

Experimental 

damage 

thresholds 

( W cm  2 ) 

Si 

1200 

1.3  x 1015 

1.10 

7.31 

x 108 * * * 

Ge 

2300 

4 x 1014 

0.79 

1.59 

x 108 

GaAs 

3500 

1 x 10 

1.40 

2.24 

CD 

O 

f-H 

X 

GaAs  P, 

x 1 — X 

(x  = .3) 

2000 

2 x 1017 

1.70 

6.25 

X 

h-* 

o 

00 

GaP 

110 

-io18 

2.20 

2.80 

9 

x 10 

1 . EXPLANATION  OF  DAMAGE  THRESHOLDS 

It  was  shown  in  Chapter  4 that  the  damage  mechanism  for 

CW  CO^  laser  was  parametric  instability  of  the  type  formulated  by 

DuBois  and  Goldman.  The  necessary  plasma  was  produced  by  photon 

induced  excitation  and  was  time  dependent.  Due  to  the  short  dura- 

“8 

tion  of  the  laser  pulse  (~3.5  x 10  sec)  in  the  present  case 

photon  induced  excitation  would  be  negligible. 

For  parametric  instability  to  grow,  an  electron  concentration 

20  -3  . 

of  1.29  x 10  cm  is  required  in  the  plasma  in  the  case  of  Si  an 
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19 

^.96  x 10  in  the  case  of  Ge.  However,  in  the  case  of  Nd : YAG 

laser  the  energy  of  a single  photon  is  1.17  eV,  which  is  higher 

than  the  band  gap  of  either  Si  (1.10  eV)  or  Ge  (0.79  eV) . It  is, 

therefore,  clear  that  band  to  band  transitions  (indirect  in  the 

case  of  Ge  and  Si)  would  be  very  dominant.  The  incident  photon 

29  3 

flux  in  the  case  of  Si  is  1.56  x 10  photons  per  second  per  cm  , 
2ft  3 

and  3.4  x 10  photons  per  second  per  cm  in  the  case  of  Ge . [The 
-2  -1 

photon  flux  in  cm  sec  has  been  divided  by  the  thickness  of  the 

sample.]  Hence,  during  the  duration  of  the  pulse,  there  will  be 

21  3 21  -3 

5.41  x 10  transitions/cm  in  Si  and  1.19  x 10  cm  in  Ge.  The 

number  of  these  transitions  is  enough  to  provide  the  required 

plasma  concentration. 

Such  band  to  band  transitions  are  energetically  forbidden  in 
GaAs,  GaAs^P^_^  (x  = 0.3)  and  GaP,  whose  band  gaps  are  higher  than 
the  photon  energy  of  1.17  eV.  But  because  of  the  high  electric 
fields  associated  with  the  incident  laser  flux,  it  is  found  that 
Zener  type  tunneling  across  the  band  gap  barrier  can  give  appre- 
ciable tunneling  rates  to  account  for  the  electron  plasma  in  GaAs, 

GaAs  P,  and  GaP. 
x 1 — x 

Zener  Tunneling  Rate 

In  this  section  an  expression  for  Zener  tunneling  [4]  rate 
using  Heisenberg  uncertainty  principal  and  WKB  approximation  will 
be  derived.  Using  this  expression,  we  shall  find  the  tunneling 


It 


certain  allowed  bands.  However,  in  the  presence  of  an  electric 

field  an  electron  can  jump  from  one  band  to  another,  provided  it 
can  travel  freely  through  a certain  distance  d,  given  by 


E 


where  E^  is  the  forbidden  band  gap  and  E is  the  applied  electric 
field. 

Quantum  mechanically,  this  implies  that  the  wave  function 
across  the  band  gap  is  exponentially  decreasing  instead  of  being 
zero  at  the  zone  boundaries.  This  is  shown  in  Figure  5-3. 


Figure  5-3.  Tunneling  through  a potential  barrier 
Applying  WKB  approximation,  the  tunneling  probability  [5] 
can  be  written  as 


...1C' 


where  g = /—  (V(x)  - E) 


and  V(x)  depends  on  the  shape  of  the  barrier  and  is  deter- 
mined both  by  the  nature  of  the  semiconductor  and  the  applied 
electric  field. 

To  calculate  the  tunneling  rate,  one  must  consider  the  natur< 

of  the  wave  function  in  the  forbidden  energy  range. 

Consider  Schroedinger  equation, 

2 

- ~ V24>  + { V(r)  — E } = 0 (5.r>) 

zm 

Try  a solution  of  the  form 


♦ -1°  „Mk-g>.r  (5.6) 

g k-g 

where  g is  a reciprocal  lattice  vector.  Substituting  Equation  5.t 
in  Equation  5.5,  produces 


- E (k ) 


\ 'V,  ► + l 

I k-g  g’ 


V a =0 
g’-g  k-g’ 


where  V(x)  has  been  expanded  in  its  fourier  series 


V(r)  = F.  V e 


i (g' -g) • r 


and  10°  = (k-g)2  (5.0) 

k-g  2m 

These  equations  are  true  for  all  values  of  g including  q 0 
and  g'  = 0.  When  the  perturbation  due  to  the  electric  field  is 
small,  then  all  ^ are  small  except  a^.  This  assumption  breaks' 
riown  near  a zone  boundary  because  of  the  fact  that  a Bragg  re- 


flection occurs.  Hence  to  construct  a stationary  wave  we  must 


;.5C- 


consider  both  the  incident  wave  and  the  diffracted  wave  on  equal 
footing. 

Ignoring  all  coefficients  except  these  two  in  Equation  5.7, 


I o 

! E - E(k) 


j <*>  + V+a->  - = 0 

1 k G k-G 


I o - j 

V + E + + " E<k>  a+  - = 0 
-G  k ( k-G  I k-G 


where 


1=1  ' T 


. _o  n2  / g \2 

and  E =2^U) 

Equations  5.10  and  5.11  will  have  a non-trivial  solution  if 
|"e°  - E(k)"l  [e°  - E 1-  |VJ2  = 0 

L *-<r  J [ k itj  1 G' 

or  h 5,2  - - E]  * IvgI2 

For  real  k,  this  equation  has  no  solution  for 


and  E 


E°  - |VG|  < E < E°  + |VQ| 


However,  this  is  the  region  of  interest  here.  Let 


k = | + < and  E = E°  + £ 


where  < and  £ are  small 


Making  this  substitution  in  Equation  5.14,  one  gets 


£2  - | V | 2 
2 = 2m  1 G1 

-2  ._o 

n 4E 


From  Equation  5.16,  we  see  that  k can  be  imaginary  if 


1 


& 


e = 


of 

Equation  5.4 

/ 2m 

(Iv0l2-.2) 

_2 

4E 

1 fi 

o 

(5.17) 


Now  using  Equation  5.17  and  Equation  5.3,  the  tunneling  proba- 
bility is  given  by 


T = exp 


- / • ■ 


dx 


(5.18) 


Now  £ = eEx  (5.19) 

where  E now  is  applied  electric  field.  From  Equation  5.2, 


l£g|2  = <IVG|2  = e2^2 


or  V„ 


2 2 2 

I 2 e E d 


(5.20) 


G1  4 

Hence  Equation  5.18  with  the  help  of  Equations  5.19  and  5.20 
bo comes 

+d/2 


T = exp  ; - 2 


’ / te  lv='2  (l  - 


-d/2 


I 


(5.21) 


Performing  the  integration  and  making  use  of  Equations  5.12 
and  5.13,  one  gets 

2 

I 

(5.22) 


T = exp 


maE 


4R2eE 


where  a is  the  lattice  constant. 

In  actual  practice  m must  be  replaced  by  the  electron  effec- 


tive mass  m* 
e 
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-■  — 


1 


' = e: 


To  obtain  the  total  tunneling  probability  one  must  multiply 
T by  the  number  of  Bragg  reflections  that  can  occur  within  a lat- 
tice period  in  one  second.  This  is  done  by  applying  Heisenberg'., 
uncertainty  principle.  If  the  electron  is  allowed  to  move  freely 
within  the  solid  through  a spacing  a before  its  momentum  is  totally 
randomized,  electric  field  has  accelerated  it  to  an  energy  of  eF.u 
so  that  the  electron  can  have  an  energy  uncertainty  of  eEa.  But 
by  the  uncertainty  principle 


AEAt  ~ h 


and  hence  the  frequency  with  which  the  Bragg  reflection  can  occur 
i s 


1 eEa 
vb  At  h 


(5.2-1) 


From  Equations  5.23  and  5.24,  the  tunneling  rate  is  given  by 


-1  eEa 

wiser  ) = -r—  exp  . - 

” I 


( m*aE 

) e g 


4R2eE 


(5.25i 


which  is  Zener's  original  formula. 


Now  tunneling  rates  for  GaAs,  GaAs^P^  ^ and  GaP  using  Equation 


..25  and  the  electric  fields  calculated  from  the  exirerimontal  duin 
age  thresholds  will  be  calculated.  The  following  expression  is 


used  to  find  the  electric  field  from  the  damage  threshold. 


The  values  for  the  effective  masses  were  taken  from  Reference 


[G] , the  lattice  constants  from  Reference  [7]  and  the  band  gaps 
from  Reference  [8] . 

The  calculated  tunneling  rates  are  listed  in  Table  5-2. 

Using  these  tunneling  rates  one  can  find  the  tunneling  probability 

— 8 

over  the  duration  of  the  pulse  which  is  3.5  x 10  seconds.  As- 
suming the  concentration  of  the  electrons  in  the  valence  band  to 

22  -3 

be  of  the  order  of  10  cm  , one  can  estimate  the  total  number  of 
electrons  per  cm^  that  can  tunnel  to  the  conduction  band  through 
this  process.  This  number  along  with  the  plasma  concentration 
required  for  the  parametric  instability  is  shown  in  Table  5-3. 

It  is  seen  from  this  table  that  the  electrons  tunneling  to  the 
conduction  band  are  enough  to  cause  the  onset  of  the  parametric 
instability. 

Calculation  of  the  Damage  Thresholds 

It  was  shown  in  Chapter  4 that  the  onset  of  parametric  in- 
stability occurs  at  a power  threshold  of 

P..  = 4irnck  T ( — \ \ [See  Equation  4.6.]  (5.27) 

th  B \up)  \ y) 

In  this  equation  damping  rate 

w 

i 

has  been  omitted,  because  over  the  duration  of  the  pulse  the  energy 
damped  due  to  this  mode  is  negligible. 

As  before,  Landau  damping  [9]  can  be  approximated  as 
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i 


1 


(5.28) 


0)  T 

P e 


“ft 


(damping  due  to  optical  phonons)  usually  lies  between  10  1 and 
-2 

10  [10].  Damping  due  to  this  mode  is  more  dominant  in  III-V 

compounds  as  compared  to  those  in  Si  and  Ge.  In  fact,  in  Si  and 
Ge  it  can  be  taken  as  1 under  the  present  conditions. 

In  Table  5-4,  the  results  of  the  calculations  of  damage 
thresholds  in  comparison  to  the  experimental  results  are  shown. 
It  is  seen  from  this  table  that  there  is  a reasonable  agreement 
between  the  theoretical  and  experimental  damage  thresholds.  The 
tact  that  the  calculated  thresholds  are  slightly  higher  than  the 
experimental  is  due  to  the  neglect  of  the  damping  mechanism 


due  to  the  acoustic  branch  of  the  lattice  spectrum. 

Shatas  et  al.[ll]  has  done  similar  calculations  of  tunneling 
rite  and  has  tried  to  explain  the  pulsed  CO^  laser  damage  in  GaAs , 
by  using  results  of  the  modified  two  stream  instability  formulated 
by  Kaw  and  Dawson  [12]  , according  to  which  the  damage  will  occur 
when 


(5.29) 


where 


eE 

E mu 
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= electron  thermal  velocity 

v^  = collision  frequency  of  the  electron 

2 2 ,22 

u>  = w + k s , 

R p e 

where  k is  the  wave  vector  of  the  interacting  laser  beam. 

Equation  5.29  can  easily  be  translated  in  terms  of  threshold 

power 
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As  pointed  out  by  Kaw  and  Dawson,  this  equation  holds  only  if 
T T\  and  this  is  not  true  in  the  present  problem  or  in  the  prob- 
lem to  which  Shatas  applied  this  result.  Shatas  et  al . incorrectly 
assumed  that  Tf  - T\  and  calculate  the  damage  thresholds  for  GaAs  for 

pulsed  CO^  laser  damage . A simple  calculation  of  damage  threshold 

9 2 

tor  GaAs  gives  9.4  x 10  W cm  and  is  two  orders  higher  than  the 
experimental  value. 

Another  type  of  photon-plasmon  instability  has  been  formulated 
by  Gersten  and  Tzoar  [13].  They  arrive  at  a threshold  condition 


_ . E lo 

c 2 4 go 
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32  it  3 W*e 


(5.31) 


Using  this  formula  for  GaAs,  one  arrives  at  a damage  threshold 
14  2 

f>t  '4  x 10  W cm  , which  is  much  higher  than  the  experimental 
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CHAPTER  6 
CONCLUSION 

Laser  induced  damage  in  semiconductors  is,  in  principle, 
similar  to  particle  and  y-ray  damage.  In  all  cases,  there  is 
lattice  damage  which  causes  degradation  in  the  electrical  prop- 
erties such  as  carrier  mobility,  carrier  concentration  and  car- 
rier lifetime.  In  spite  of  this  apparent  similarity,  there  are 
some  very  important  differences  that  must  be  noted.  While  the 
amount  of  damage  is  directly  proportional  to  the  incident  par- 
ticle flux  or  the  y-ray  intensity,  a sharply  defined  power  level 
exists  at  which  laser  damage  can  occur.  Also,  particle  and  y-ray 
damage  anneals  faster  than  laser  induced  damage.  This  proves  that 
the  fundamental  damage  mechanism  in  laser  damage  is  different 
than  in  the  particle  and  the  y-ray  damage. 

The  experimental  damage  thresholds  are  found  to  be  different 

at  different  laser  wavelengths.  The  various  semiconductors,  e.g. 

Si,  Ge , GaAs , GaAs  P,  and  GaP , when  radiatdd  with  CW  C0_  laser 
x l-x  2 

2 

radiation  were  all  damaged  at  power  thresholds  of  a few  kW/cm*' , 
the  irradiation  time  being  0.10  second.  When  a systematic  analysis 
of  different  damage  mechanisms  is  made,  it  is  found  that  the  band 
to  band  transitions,  the  multiphoton  absorption,  the  free  carrier 
absorption  or  the  lucky  electron  theory  of  avalanche  cannot  ex- 
plain the  experimentally  observed  damage  thresholds.  This  leaves 
one  with  the  plasma  instabilities.  It  is  found  that  semiconductors 
ran  be  treated  as  containers  of  electron  plasma  that  behave  simi- 
larly to  gaseous  plasmas.  The  observed  damage  thresholds  can  be 
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satisfactorily  explained  by  applying  the  theory  of  parametric  in- 
stability of  the  type  formulated  by  DuBois  and  Goodman.  This  in- 
stability  causes  anomalous  increases  in  the  semiconductor  resis- 
tivity causing  anomalous  absorption  as  suggested  by  Kaw  and  Dawson. 

The  pre-condition  for  the  onset  of  the  parametric  instability 

17  18 

at  10.6  pm  is  a plasma  concentration  of  10  to  10  electrons 
-3  2 

cm  . Also  A (a  parameter  defined  in  Chapter  2)  should  be  less 
than  1.0.  The  latter  condition  is  found  to  be  easily  satisfied, 
while  the  required  plasma  concentration  is  achieved  by  photon  in- 
duced excitation  of  valence  electrons  to  the  conduction  band  through 

a continuous  distribution  of  Shockley  surface  states.  The  time  re- 

17  18 

quired  for  the  plasma  concentration  to  build  up  to  10  to  10 
cm  ^ is  of  the  order  of  a few  milliseconds.  There  is  enough  time 
available  for  this  to  happen  in  the  case  of  CW  C0„  radiation  dam- 
age.  This  mechanism  cannot  account  for  a solid  state  plasma  con- 
centration of  10^  to  10^®  electrons  cm  ^ required  for  the  onset 
of  parametric  instability  at  1.06  pm  because  of  the  very  short 
duration  of  the  pulse.  The  experimentally  determined  surface  dam- 
age thresholds  at  1.06  pm  for  nanosecond  pulses  were  of  the  order 
8 9 2 

of  10  to  10  W/cm  . These  damage  thresholds  were  higher  for 

higher  band  gap  semiconductors.  As  the  energy  of  a single  photon 

at  1.06  pm  is  1.17  eV  and  the  incident  photon  flux  is  of  the  order 
20  21  -3 

of  10  to  10  cm  , the  required  plasma  concentration  can  be 
achieved  through  direct  or  indirect  band  to  band  transitions  in 
the  cases  of  Ge  and  Si.  Such  transitions  are  energetically  for- 
biden  in  GaAs , GaAs^P ^ and  GaP  because  of  the  higher  band  gaps 
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of  these  semiconductors.  A simple  calculation  shows  that  the 

9 -2 

longitudinal  electric  field  associated  with  a flux  of  10  W cm 

5 

is  of  the  order  of  10  V/cm.  The  Zener  tunneling  rates  at  this 

value  of  the  electric  field  are  found  to  be  of  the  order  of  10^ 

6 “ 1 »2  — 3 

to  10  sec  , giving  a tunneling  probability  of  10  to  10  over 

22 

the  duration  of  the  pulse.  Because  there  are  about  10  electrons 

cm  3 in  the  valence  band  of  any  one  of  these  semiconductors,  the 

19 

number  of*electrons  that  tunnel  to  the  conduction  band  is  10  to 

20  -3 

10  cm  , thus  giving  the  required  plasma  concentration  for  the 
onset  of  the  parametric  instability  at  1.06  ym. 

In  practice damage  is  found  to  occur  both  at  the  surface  and 
within  the  bulk  of  the  semiconductor.  This  creates  defects  in  the 
form  of  interstitial  and  vacancy  pairs,  giving  rise  to  the  local- 
ized traps  within  the  band  gap.  These  traps  are  found  to  be  more 
active  at  lower  temperatures  than  at  room  temperature,  as  seen  from 
the  experimental  Hall  mobility  versus  temperature  curves.  Due  to 
this  reason,  these  defects  can  be  considered  as  scattering  centers 
similar  to  ionized  impurities.  Both  in  Germanium  and  Silicon,  it 
is  found  that  the  mobility  decrease  due  to  CW  CO^  laser  damage  is 
about  30  to  50  percent,  and  also  there  is  significant  carrier  re- 
moval in  these  semiconductors.  In  this  respect,  therefore,  the 
laser  damage  is  similar  to  the  particle  and  the  y-ray  damage. 

These  results  on  mobility  and  carrier  concentration  can  be  ex- 
plained by  applying  the  model  of  James  and  Lark-Horovitz,  assuminq 
the  creation  of  localized  levels  at  0.18  eV  below  the  conduction 
band  in  the  case  of  Ge  and  0.36  eV  in  the  case  of  Si. 
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